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ABSTRACT 


Detailed  experimental  investigation  of  second  generation  controlled-diffusion 
compressor  stator  blades  at  design  inlet-flow  angle  was  performed  in  a  low-speed  cascade 
wind  tunnel  using  various  experimental  methods.  Surface  pressure  measurements  were 
obtained  using  three  instrumented  blades,  from  which  coefficients  of  pressure  were 
calculated.  Laser-Doppler  velocimetry  was  used  to  characterize  the  flow  in  the  inlet,  in 
the  passage  between  two  blades,  in  the  boundary  layer  of  the  blades,  and  in  the  wake.  A 
five-hole  pressure  probe  was  used  to  determine  the  loss  coefficient  and  the  axial-velocity- 
density  ratio  of  the  flow  through  the  cascade.  Although  the  blades  produced  significant 
lift,  separated  flow  was  discovered  on  the  suction  side  of  the  blades  at  approximately  fifty 
percent  axial  chord,  which  showed  that  the  design  was  not  totally  successful.  All  the 
experimental  measurements  were  performed  at  an  inlet  flow  Mach  number  of  0.22  and  a 
Reynolds  number,  based  on  chord  length,  of  640,000. 

Experimental  blade-surface  pressure  coefficients  were  compared  with  values 
predicted  using  a  computational  fluid  dynamics  code.  These  initial  predictions  did  not 
match  well  with  the  experimental  results. 
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1.  INTRODUCTION 


A.  BACKGROUND 

Increased  performance  requirements  for  aircraft  gas  turbine  engines  necessitate 
continuous  progress  in  gas  turbine  engine  research  and  development.  Engine  compressor 
stall  and  off-design  behavior  continue  to  limit  the  performance  of  military  aircraft  gas 
turbine  engines.  Compressor  stall  can  lead  to  degradation  in  engine  performance  or  even 
total  loss  of  engine  power,  which  could  result  in  mission  failure  or  aircraft  loss. 
Additionally,  the  increased  thrust-to-weight  ratio  engines  needed  for  present  and  future 
military  and  civil  aircraft  require  improved  compressor  design  for  increased  performance 
at  the  same  or  reduced  weight. 

Compressor  blading  has  progressed  over  the  years  from  the  use  of  NACA-65 
series.  Double  Circular  Arc  (DCA)  and  Multiple  Circular  Arc  (MCA)  bladeshapes  to  the 
design  of  Controlled-Diffusion  (CD)  blading.  CD  blading  was  developed  to  control 
diffusion  on  the  suction  side  of  the  airfoil,  thus  avoiding  boundary-layer  separation  at  the 
engine  design  point.  Transonic  design  allowed  shock-free  operation  in  the  transonic  range 
which  increased  compressor  and  engine  efficiency.  The  resulting  increase  in  efficiency 
could  be  used  to  improve  engine  performance,  or  the  same  engine  performance  could  be 
achieved  using  fewer  compressor  blades,  allowing  a  reduction  in  engine  weight. 

The  CD  compressor  blades  investigated  in  the  present  study  were  designed  by 
Thomas  F.  Gelder  of  NASA  Lewis  Research  Center  (LeRC)  [Ref  1].  The  compressor 
stator  profiles  were  Stator  67B  blades,  which  together  with  Rotor  67  comprised 
Compressor  Stage  67B.  Compressor  Stage  67B  was  previously  studied  experimentally  at 
the  NASA  LeRC  compressor  test  facility  [Ref  1].  The  Stator  67B  blades  were  second 
generation  CD  blades  which  were  designed  as  an  improvement  to  the  Stator  67A  CD 
blades.  The  Stator  67A  blades,  designed  by  Nelson  Sanger  [Ref.  2] ,  coupled  with  Rotor 
67  to  form  Stage  67 A.  Ten  midspan  Stator  67B  blade  profiles  were  machined  from 
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aluminum  and  installed  in  the  Naval  Postgraduate  School  (NPS)  Turbopropulsion 
Laboratory  (TPL)  Low-Speed  Cascade  Wind  Tunnel  (LSCWT)  for  the  present  study. 
Pneumatic  probe  measurements  were  made  upstream  and  downstream  of  the  blading. 
Laser-Doppler  velocimetry  (LDV)  measurements  were  performed  upstream  of  the  blade 
row,  in  the  passage  between  the  blades,  in  the  boundary  layers,  and  downstream  in  the 
wake,  at  the  design  inlet-flow  angle.  A  quasi-three  dimensional  computational  fluid 
dynamics  (CFD)  code.  Rotor  Viscous  Code  Quasi-3-D  (RVCQ3D)  Version  300,  was  also 
used  for  comparison  with  the  experimental  results. 

The  Stator  67A  compressor  blading  was  the  focus  of  a  decade  of  research  at  the 
NPS  TPL.  Koyuncu  [Ref  3]  studied  the  performance  at  various  incidence  angles  using  a 
five-hole  probe.  Dreon  [Ref  4]  improved  the  accuracy  in  loss  and  pressure  measurement 
near  design  conditions.  Elazar  [Ref  5]  performed  LDV  measurements  in  the  boundaiy 
layer,  in  the  blade  passage  and  in  the  wake  at  three  incidence  angles.  Murray  [Ref  6] 
upgraded  the  LDV  system  and  measured  the  flow  in  the  wake  near  stall.  Classick  [Ref  7] 
upgraded  the  pressure  probe  data  acquisition  software  and  performed  surveys  near  stall. 
Armstrong  [Ref  8]  further  improved  the  data  acquisition  system  and  compared  flow  at 
high  incidence  angles  over  a  blade  with  a  modified  leading  edge  and  an  unmodified  blade. 
Hobson  and  Shreeve  [Ref  9]  used  LDV  to  study  the  flow  at  very  high  incidence.  Ganiam 
[Ref  10]  performed  LDV  measurements  at  stall.  Williams  [Ref  11]  verified  the  LDV 
measurements  at  stall  and  obtained  a  prediction  using  computational  fluid  dynamics 
(CFD). 

B.  PURPOSE 

The  objective  of  the  present  study  was  to  install  the  Stator  67B  midspan  blade 
sections  in  the  NPS  cascade  and  to  examine  the  flow  through  the  blading  at  the  design 
inlet-flow  angle  using  LDV  and  pressure  probe  measurements.  The  intent  was  first  to 
obtain  a  thorough  understanding  of  the  flow  at  the  design  inlet-flow  angle,  and  to  validate 
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the  design  itself.  Future  studies  of  Stator  67B  cascade  blading  at  higher  inlet-flow  angles 
will  determine  the  off-design  and  stalling  behavior. 
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II.  TEST  FACILITY  AND  INSTRUMENTATION 


A.  LOW-SPEED  CASCADE  WIND  TUNNEL 

The  NPS  Turbopropulsion  Laboratory  Low-Speed  Cascade  Wind  Tunnel  Building 
is  shown  in  Figure  1.  Elazar  [Ref.  5]  thoroughly  documented  the  uniformity  of  tunnel 
flow  conditions  and  the  periodicity  of  the  flow  in  the  cascade  test  section  with  20  Stator 
67A  blades  in  the  cascade  at  40  (design),  43  and  46  degrees  inlet-flow  angle. 


B.  TEST  SECTION 

A  schematic  of  the  cascade  is  shown  in  Figure  2.  Ten  CD  blades  with  elliptical 
leading-  and  trailing-edges  were  scaled  from  the  midspan  section  of  Stator  67B  [Ref  1] 
and  machined  from  aluminum  using  numerically-controlled  machining  methods.  Figure  3 
shows  the  profile  of  the  blade.  Table  1  contains  the  machine  coordinates  used  to 
manufacture  the  blades.  Three  blades  were  instrumented  with  pressure  taps  which  were 
machined  into  fine  metal  tubing  laid  into  the  blade  surface.  A  fiilly-instrumented  blade 
containing  42  pressure  taps  was  installed  at  blade  location  number  6  [Figure  2].  Two 
blades  numbered  2  and  8  were  partially  instrumented  with  eight  pressure  taps  each.  Figure 
4  shows  the  pressure  tap  locations  on  the  instrumented  blades.  Additionally,  blades  3  and 
4  were  treated  with  a  black  anodized  coating  to  minimize  laser  light  scatter  for  the  LDV 
measurements.  The  LDV  test  passage  between  blades  3  and  4  is  shown  in  Figure  5  with 
the  locations  of  the  13  measurement  stations  shown  as  fractions  of  the  axial  chord  (Cac). 
Table  2  contains  a  summary  listing  of  the  geometrical  parameters  of  the  cascade  test 
section. 

The  ten  blades  were  installed  in  the  test  section  using  brass  shims  for  alignment 
perpendicular  to  the  cascade  endwalls.  A  digital  inclinometer  with  an  accuracy  of  0.1 
degrees  was  used  to  set  the  stagger  angle  of  the  blades.  The  tunnel  inlet  sidewalls  were 
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adjusted  to  the  design  inlet-flow  angle,  Pi  =  36.3  degrees.  The  inlet  guide  vanes  were 
adjusted  until  the  mean  inlet  flow  to  the  test  section  at  Station  1  was  equivalent  to  the 
design  inlet-flow  angle.  The  final  tunnel  adjustment  was  that  of  the  tail  board  setting 
angles  which  were  adjusted  to  give  uniform  wall  static  pressures  downstream  of  the 
blading  [Figure  2]. 


C.  INSTRUMENTATION 

1.  Pneumatic  Data  Acquisition  System 

Two  different  pneumatic  data  acquisition  systems  were  used  with  a  traversing  five- 
hole  probe  to  determine  the  loss  coefficient  (co),  the  axial-velocity-density  ratio  (AVDR), 
and  blade  surface  pressure  coefficients  (Cp).  The  first  system  was  a  computer-controlled 
automated  data  acquisition  system  and  software  documented  by  Classick  [Ref  7]  and 
modified  by  Armstrong  [Ref  8].  The  Hewlett-Packard  data  acquisition  system  hardware 
and  software  “ACQUIRE”  and  “LOSS”  are  fully  described  by  Armstrong  [Ref  8]. 

The  second  data  acquisition  system  consisted  of  a  Scanivalve  Digital  Interface  Unit 
(SDIU  Mk5),  a  Scanivalve  Controller  (CTLR2/S2-S6),  a  Hewlett-Packard  digital 
voltmeter  (HP  3437A)  and  a  48  channel  Scanivalve  rotary  pressure  scanner.  Loss  and 
AVDR  calculations  were  performed  using  a  personal  computer  with  the  spreadsheet 
software  Microsoft  “Excel”.  The  five-hole  probe  and  the  five-hole  probe  traverse  used  for 
both  systems  are  shown  in  the  photographs  in  Figures  6  and  7.  Figure  8  shows  a 
schematic  of  the  second  data  acquisition  system.  Plenum  pressure  and  temperature 
sensors  were  the  same  as  those  used  in  the  first  system.  A  Prandtl  probe  located  upstream 
of  the  test  section  was  used  as  the  total  pressure  reference  for  the  second  system. 

Prior  to  implementation  of  either  of  the  two  acquisition  systems,  a  series  of  surface 
pressure  measurements  was  performed  using  banks  of  water  manometers.  The  partially- 
instrumented  blades,  mounted  at  blade  position  numbers  2  and  8  in  the  test  section,  were 
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connected  to  manometers  for  an  initial  observation  of  the  pressure  distribution  over  the 
blades. 


2.  Laser  Doppler  Velocimeter 

A  four  beam,  two  color  TSI  Model  9100-7  LDV  system  was  used  for  all  surveys. 
Elazar  [Ref  5]  thoroughly  described  the  LDV  system,  including  the  laser  model,  optics, 
atomizer  and  seeding,  and  data  acquisition.  Murray  [Ref  6]  described  the  fully  automated 
traverse  mechanism.  The  LDV  and  traverse  systems  were  controlled  by  a  personal 
computer.  Experimental  data  were  processed  using  TSI  FIND  software  on  the  personal 
computer.  The  LDV  and  traverse  systems  are  shown  in  Figure  9. 
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Figure  2.  Low-Speed  Cascade  Wind  Tunnel. 
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Figure  5.  Survey  Station  Numbering  and  Position 


Blade  Type 

Stator  67B  Controlled-Difflision 

Number  of  Blades 

10 

Blade  Spacing 

152.4  mm 

Chord 

127.14  mm 

Solidity 

0.834 

Thickness/Chord 

0.05 

Setting  Angle 

16.3°  ±  0.1° 

Span 

254.0  mm 

Table  2.  Test  Section  Data. 
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Port  #1  (P  Atmosperic) 
Port  #2  (P  Calibration) 

Port  #3  (P  Plenum) 

Port  m  (P  Wall  Static) 
Port  #5  (PI  5-Hole  Probe) 
Port  #6  (P2  5-Hole  Probe) 
Port  #7  (P3  5-hole  Probe) 
Port  #8  (P4  5-Hole  Probe) 
Port  #9  (P5  5-Hole  Probe) 
Port  #1 0  (P  Prandtl  Total) 


HP  3437A  Digital  Port  #1 1  (P  Prandtl  Static) 

VoUmeter 


Scanivalve 

Controller 

CTLR2/S2'S6 


Figure  8.  Manual  Five-Hole  Probe  System. 
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III.  EXPERIMENTAL  PROCEDURE 


A.  SURFACE  PRESSURE  MEASUREMENTS 

The  tunnel  was  started  and  allowed  to  reach  an  equilibrium  operating  temperature 
prior  to  all  surveys.  Once  the  tunnel  reached  operating  temperature,  the  plenum  chamber 
pressure  was  set  and  maintained  at  12  inches  of  water.  Twelve  inches  of  water  for  the 
plenum  pressure  gave  an  average  inlet  Mach  number  of  0.22  and  an  average  Reynolds 
number  based  on  chord,  Rcc,  of  approximately  640,000. 

1.  Water  Manometer  System 

Prior  to  the  automated  data  acquisition  system  experiments,  surface  pressure 
measurements  were  obtained  using  banks  of  water  manometers  connected  to  the  two 
partially-instrumented  blades.  The  water  manometers  were  manually  recorded  to  an 
accuracy  of  0.1  inches  of  water.  Plenum  total  pressure  and  temperature,  wall  static 
pressure,  Prandtl  probe  total  and  static  pressures,  and  local  atmospheric  pressure  were 
recorded  in  addition  to  the  blade  pressure  taps.  Calculations  of  the  coefficients  of 
pressure,  Cp,  based  on  the  Prandtl  probe  total  and  static  pressures,  were  performed  using 
the  software  Microsoft  “Excel”. 

2.  HP  Automated  Data  Acquisition  System 

The  HP  data  acquisition  system  was  initialized,  and  the  Scanivalve  was  calibrated 
using  a  water  manometer  prior  to  each  experimental  run.  A  U-tube  water  manometer  and 
a  regulated  compressed-air  supply  were  used  as  the  calibration  reference  for  both  of  the 
Scanivalve  rotary  pressure  scanners.  The  five-hole  probe  was  positioned  away  from  the 
three  instrumented  blades  to  prevent  wake  interference.  The  program  “ACQUIRE” 
sampled  and  recorded  pressures  using  a  Scanivalve  rotary  pressure  scanner,  and  then 
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calculated  Cp  for  each  blade  pressure  tap  based  on  the  plenum  total  pressure,  as  described 
by  Armstrong  [Ref  8], 


B.  FIVE-HOLE  PRESSURE  PROBE  MEASUREMENTS 

The  five-hole  probe  traverse  mechanism  was  mounted  upstream  of  the  test  section 
at  the  tunnel  lower  survey  location  for  the  loss  and  AVDR  measurements.  The  five-hole 
probe  was  then  mounted  on  the  traverse  mechanism  and  centered  at  the  blade  midspan 
position  in  the  tunnel.  A  U-tube  water  manometer  was  used  with  a  regulated  compressed- 
air  supply  to  calibrate  the  Scanivalve  rotary  pressure  scanner  prior  to  making  the 
measurements.  The  five-hole  probe  yaw  transducer  was  also  calibrated  using  a  digital 
inclinometer.  The  five-hole  probe  was  positioned  in  the  pitchwise  direction  and  aligned 
with  the  tunnel  flow  by  balancing  the  probe  in  yaw  using  a  U-tube  water  manometer  prior 
to  taking  pressure  data.  The  five-hole  probe  calibration  coefficients  developed  by 
Armstrong  [Ref  8]  were  used  by  the  data  analysis. 

1.  HP  Automated  Acquisition  System 

The  HP  automated  acquisition  system  was  used  to  determine  loss  and  AVDR  for 
the  blade  section  as  outlined  by  Armstrong  [Ref  8].  The  five-hole  probe  was  positioned 
in  both  the  lower  and  the  upper  traverse  positions.  Data  were  taken  over  a  full  passage 
width  of  154  mm  in  increments  of  2.54  mm,  5.08  mm  and  12.77  mm,  respectively,  for  the 
surveys  using  Armstrong’s  “ACQUIRE”  program  [Ref  8].  The  “LOSS”  program  was 
then  used  to  calculate  loss  and  AVDR  for  each  combination  of  upstream  and  downstream 
data,  as  described  by  Armstrong  [Ref  8]. 

2.  Manual  Loss  Coefficient  and  AVDR  Measurements 

The  SDIU,  Scanivalve  controller,  Scanivalve  rotary  pressure  scanner  and  HP 
digital  voltmeter  were  used  with  the  five-hole  probe  to  determine  the  loss  coefficient  and 
AVDR.  The  SDIU  and  controller  were  used  to  step  the  Scanivalve  to  each  port  for 
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measurements,  and  the  pressure  data  were  read  from  the  digital  voltmeter  and  recorded 
into  a  spreadsheet  on  a  personal  computer  for  data  reduction.  Plenum  total  pressure  and 
temperature,  atmospheric  pressure,  wall  static  pressure,  and  Prandtl  probe  total  and  static 
pressures  were  recorded.  Probe  yaw  angle  was  also  recorded.  The  spacing  interval  for 
the  survey  was  2.54  mm  over  a  total  survey  width  of  152.4  mm  with  a  total  of  61  data 
positions  being  recorded. 

Appendix  A  contains  the  formulas  which  were  used  for  the  loss  calculations.  The 
calculation  method  was  as  follows.  First,  the  dimensionless  velocity  coefficients,  3,  were 
determined  from  the  five-hole  probe  pressure  data.  Next,  the  dimensionless  velocities,  X, 
were  determined  using  the  polynomial  coefficients  for  the  five-hole  probe.  The 
dimensionless  reference  velocity,  Xref,  was  determined  using  the  pressures  from  the  Prandtl 
probe.  The  reference  flow  (massflux)  functions,  K,  were  calculated  for  both  the  upstream 
and  the  downstream  data.  AVDR  was  determined  using  a  Romberg  numerical  integration 
routine  from  the  values  of  K.  The  loss  coefficient,  co,  was  determined  from  values  of  CpK, 
using  the  same  type  of  numerical  integration  scheme. 


C.  LDV  MEASUREMENTS 

1.  LDV  Probe  Volume  Alignment  and  Reference 

The  need  for  a  precise  reference  for  the  LDV  probe  volume  necessitated  the  design 
of  an  alignment  tool.  A  photograph  of  the  tool,  which  was  manufactured  out  of 
aluminum,  is  shown  in  Figure  10,  and  the  tool  machine  drawing  with  dimensions  is 
presented  in  Appendix  B.  The  tool  was  positioned  on  blades  3  and  4  in  the  spanwise 
direction  by  placing  the  two  end  pieces  against  the  back  of  the  test  section  wall.  The  two 
endpieces  secured  the  tool  tightly  against  the  trailing  edge  of  both  blades  using  the 
machine  screw  on  the  right-hand  side,  which  ensured  consistent  positioning.  Three  laser 
alignment  holes  were  machined  in  the  tool’s  center  block,  each  with  a  0.3302  mm 
diameter.  The  laser  probe  volume  was  focused  through  one  of  these  small  holes,  and  the 
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automated  traverse  table  (TSI  Model  9500)  was  initialized  to  the  coordinates  of  the  hole. 
This  alignment  process  was  performed  for  each  survey. 

The  zero  reference  point  was  determined  using  the  blade  geometry  at  the  design 
setting  angle.  The  reference  point  was  defined  by  the  normal  intersection  of  two  lines 
which  were  tangent  to  the  minimum  x  and  y  points  of  the  leading-edge  ellipse  of  blade 
number  3.  At  this  intersection,  x  and  y  were  defined  to  be  zero.  The  spanwise  mid-point 
of  the  blade,  at  127  mm,  was  used  as  the  third  zero  reference  point,  and  all  measurements 
were  performed  at  midspan. 

2.  Inlet  Guide  Vanes  Adjustment 

Initial  LDV  inlet  surveys  were  performed  over  a  254  mm  distance,  or  nearly  two 
blade  passage  widths,  at  Station  1  to  allow  proper  adjustment  of  the  inlet  guide  vanes. 
The  mean  inlet  flow  angle  was  computed  from  survey  data  and  the  inlet  guide  vanes  were 

adjusted  until  this  angle  converged  to  within  0.1  degree  of  the  design  angle  of  36.3 

degrees. 

3.  LDV  Surveys 

A  total  of  28  surveys  were  completed,  including  at  least  one  survey  of  each  of  the 
thirteen  stations  and  four  boundary  layer  surveys.  All  surveys  were  taken  by  collecting 
1000  data  points  at  each  location  across  the  station,  with  the  exception  of  the  Station  1 
surveys,  where  3000  data  points  were  taken.  Figure  1 1  shows  the  location  of  the  station 
surveys  which  will  be  discussed  in  the  body  of  this  report.  A  listing  of  all  experimental 
runs  performed  is  included  in  Appendix  C. 

The  laser  optics  were  configured  identically  throughout  this  study,  with  the  514.5 
nm  green  beam  measuring  the  axial  (or  vertical)  velocity  component,  U,  and  with  the  488 
nm  blue  beam  measuring  the  tangential  (or  horizontal)  velocity  component,  V.  Frequency 
shifting  of  5  MHz  was  performed,  as  outlined  by  Ganaim  [Ref  10],  to  detect  reverse  flow 
velocities.  Alignment  was  performed  prior  to  each  measurement  using  the  laser  alignment 
tool,  as  described  above. 
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Ambient  pressure,  and  plenum  total  pressure  and  plenum  total  temperature  were 
recorded  for  each  survey.  The  velocity  at  the  inlet  to  the  test  section  was  used  as  the 
tunnel  reference  velocity,  Vref.  The  reference  velocity  was  used  to  nondimensionalize  the 
velocity  measurements.  A  series  of  calibration  surveys  at  six  plenum  pressures  and 
temperatures  were  performed  and  a  least-squares  curve  fit  was  applied  to  the  data  to 
determine  a  calibration  curve  for  Vref.  A  Newton-method  iteration  algorithm  was  used  to 
determine  Vref  for  each  survey  using  plenum  total  pressure,  plenum  total  temperature  and 
atmospheric  pressure.  A  FORTRAN  code  “C ALIBI”  was  used  to  calculate  Vref,  and  the 
FORTRAN  code  listing  and  input  and  output  data  files  are  contained  in  Appendix  D. 

TSI  Flow  Information  Display  (FIND)  Version  4.0  software  was  used  to  acquire 
and  analyze  all  LDV  data.  Velocity,  turbulence  intensity,  Reynolds  stress  and  Reynolds 
stress  correlation  coefficient  information  were  processed  using  FIND,  and  then  further 
non-dimensionalized  using  the  inlet  flow  reference  velocity,  Vref,  to  allow  comparison  with 
data  taken  at  different  LDV  configurations  (i.e.  LDV  axis  not  perpendicular  to  the  optical 
access  window). 

a.  Inlet  Surveys 

The  inlet  flow  region  was  surveyed  across  Stations  1,  2  and  3,  in  the  far- 
and  near-upstream  regions.  The  laser  system  was  positioned  horizontally  for  the  Station  1 
and  2  measurements.  Station  3  measurements  were  conducted  first  with  no  pitch,  for  a 
region  comprising  93  percent  of  passage  width,  and  then  the  laser  was  pitched  upward 
four  degrees,  allowing  the  full  passage  to  be  studied.  The  potential  problems  due  to 
pitching  and  yawing  of  the  LDV  system  were  discussed  by  Hobson  and  Shreeve  [Ref  9]. 
Their  analysis  showed  a  maximum  spatial  error  fi-om  probe  volume  orientation  to  be  0.3 
mm,  which  is  the  measurement  volume  minimum  diameter. 

b.  Passage  Surveys 

Passage  surveys  were  conducted  at  Stations  4  through  10.  The  laser  was 
horizontal  for  all  passage  surveys.  One  survey  taken  at  Station  7  was  performed  with  a 


23 


laser  yaw  angle  of  4  degrees,  but  no  pitch.  Initial  and  final  measurement  boundaries  were 
determined  for  the  passage  surveys  by  experimenting  with  the  laser  position  until  an 
adequate  data  rate  was  achieved.  Interference  degraded  the  laser  beams  when  positioned 
too  close  to  the  blade  surface,  thus  diminishing  the  data  rate  of  the  LDV  system  to 
unusably  low  levels.  The  boundaries  of  the  surveys  were  set  to  exclude  regions  of 
interference. 


c.  Boundary  Layer  Surveys 

Boundary  layer  surveys  were  performed  at  one  station  on  the  pressure  side 
and  three  stations  on  the  suction  side  of  the  blade.  The  boundary  layer  surveys  on  the 
suction  surface  of  blade  3  were  performed  at  Stations  5,  7  and  9.  The  boundary  layer 
survey  on  the  pressure  surface  of  blade  4  was  performed  at  Station  8.  Prior  to  the 
experiment,  the  laser  alignment  tool  was  used  to  align  the  laser  at  the  measurement  pitch 
and  yaw  angles.  The  laser  pitch  and  yaw  angles  used  for  the  surveys  were  determined 
experimentally  by  finding  the  angle  which  allowed  the  closest  laser  probe  volume 
positioning  to  the  blade  surface  while  giving  an  adequate  data  rate.  The  pitch  and  yaw 
angles  which  were  used  for  the  boundary  layer  surveys  can  be  found  in  the  summary  in 
Appendix  C.  The  boundary  layer  surveys  were  performed  along  a  line  normal  to  the  blade 
surface  at  that  station. 

d.  Wake  Surveys 

Wake  surveys  of  various  detail  were  performed  at  Stations  11,  12,  and  13. 
Station  1 1,  the  near  wake,  and  Station  13,  the  far  wake,  are  presented  in  the  next  section 
of  this  report.  The  laser  pitch  angle  was  5  degrees  down  for  the  Station  1 1  wake  surveys. 
The  laser  was  horizontal  for  the  Station  13  surveys.  The  incremental  distance  between 
measurement  positions  was  reduced  over  a  series  of  three  surveys  at  each  of  the  two 
stations,  from  a  full  passage  width  and  large  increments  to  a  smaller  passage  width  and 
smaller  increments.  Again,  the  laser  was  aligned  prior  to  each  survey  using  the  alignment 
tool. 
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D.  COMPUTATIONAL  FLUID  DYNAMICS 


A  quasi-three-dimensional  computational  fluid  dynamics  program  was  used  to 
analyze  the  flow  through  the  cascade.  The  FORTRAN  code.  Rotor  Viscous  Code  Quasi- 
3-D  (RVCQ3D)  Version  300,  was  written  by  Rodrick  V.  Chima  of  NASA  LeRC  [Ref 
12],  A  340  X  49  C-type  grid  was  generated  using  the  FORTRAN  program  “GRAPE” 
[Ref  13]  based  on  blade  manufacturing  dimensions  [Table  1].  Appendix  E  contains  the 
input  data  used  for  “GRAPE”.  The  computational  grid  generated  by  “GRAPE”  is  shown 
in  Figure  12.  RVCQ3D  had  the  option  to  use  three  different  turbulence  models:  the 
Baldwin-Lomax  model,  the  Cebeci-Smith  model,  and  the  Wilcox  k-omega  model;  all  were 
used.  The  design  inlet-flow  angle,  Mach  number,  Reynolds  number  and  pressure  ratio 
were  other  inputs  required  to  run  the  program.  The  code  was  run  for  various  pressure 
ratios  until  the  design  inlet  flow  angle  of  36.3  degrees  and  the  Mach  number  of  0.22  were 
met.  [See  Appendix  E  for  a  sample  of  the  input  data  used  for  RVCQ3D].  A  FORTRAN 
program  modified  by  Williams  [Ref  11]  called  “PCP”  was  used  to  extract  blade  surface 
pressure  distribution  data  from  the  solution  file  in  order  to  plot  coefficients  of  pressure, 
Cp.  The  program  “PCP”  is  also  included  in  Appendix  E. 
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IV.  RESULTS  AND  DISCUSSION 


A.  REFERENCE  CONDITIONS 

The  results  of  the  experimental  data  will  be  discussed  in  the  following  sections. 
Nominally,  all  surveys  were  performed  at  a  Mach  number  of  0.22  and  a  Reynolds  number 
based  on  blade  chord  of  640,000. 


B.  BLADE  SURFACE  PRESSURE  DISTRIBUTION 

Figure  13  shows  the  results  of  the  blade  surface  pressure  distribution 
measurements  in  terms  of  the  coefficient  of  pressure,  Cp,  plotted  along  the  blade  chord  at 
various  positions  given  by  the  ratio  ^/c.  Blades  2  and  8  were  partially-instrumented  with 
eight  pressure  taps  each.  The  subscript  “m”  in  the  figure’s  legend  denotes  the  manometer 
readings  used  on  blades  2  and  8.  Blade  6  contained  42  pressure  taps,  including  a  leading- 
edge  tap  and  a  trailing-edge  tap  and  20  taps  on  each  of  the  pressure  and  suction  surfaces. 
Two  taps  were  unusable  on  the  pressure  surface  of  blade  6,  so  a  total  of  40  pressure 
measurements  was  made.  Gelder’s  design  intent  for  Stator  67B  [Ref.  1]  is  shown  plotted 
with  the  solid  black  line  in  Figure  13. 

The  results  from  the  different  blades  and  different  data  acquisition  methods 
showed  good  agreement.  The  Cp  calculations  from  both  the  water  manometer  readings 
and  the  data  acquisition  system  produced  very  similar  results.  The  Cp  results  for  the  blade 
pressure  surface  were  very  tightly  grouped  together  for  all  three  blades  and  both 
measurement  techniques.  The  blade  suction-side  results  showed  more  loosely  grouped 
data  for  blade  numbers  2  and  8  for  the  two  data  acquisition  methods.  Blade  number  2  had 
an  unusable  pressure  tap  at  0.95  ^/c,  and  also  showed  a  partially  clogged  pressure  tap  at 
0.69  ^/c.  The  suction  side  data  from  blades  2  and  8  corresponded  well  to  the  data  on 
blade  6  at  all  other  locations. 
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The  leading-edge  pressure  tap  on  blade  number  6  showed  the  stagnation  point  to 
be  very  near  the  leading-edge,  since  the  leading-edge  value  was  Cp  =  0.92.  Stagnation 
would  give  unity.  This  was  a  confirmation  of  the  design  inlet-flow  angle  which  resulted  in 
zero  incidence  at  the  leading-edge.  The  trailing-edge  pressure  tap  showed  the  value  of  Cp 
to  be  0. 10,  or  very  close  to  zero. 

Gelder’s  experimental  results  for  the  Stage  67B  [Ref  1]  were  similar  to  the 
present  results  in  terms  of  the  relation  of  the  experimental  data  to  the  design  intent.  On 
the  pressure  side  of  the  blade,  Gelder’s  experimental  Cp  data  were  lower  in  magnitude 
than  the  design  value,  similar  to  the  NFS  results.  On  the  suction  side  of  the  blade,  the 
Gelder  experimental  data  had  smaller  (negative)  magnitudes  than  design  until 
approximately  0.2  ^/c,  at  which  point  the  experimental  data  had  larger  magnitudes  than 
the  design  data. 

The  experimental  Cp  profile  seemed  to  show  that  the  design  goal  of  controlled 
diffusion  was  achieved.  On  the  suction  side  of  the  blade,  Cp  decreased  (negatively)  very 
sharply  from  the  leading-edge  of  the  blade  as  the  flow  accelerated.  Cp  then  decreased 
more  gradually  until  at  nearly  midchord,  where  Cp  began  a  near-linear  increase,  to  arrive 
close  to  zero  at  the  trailing-edge.  The  near-linear  increase  showed  the  diffusion  or 
deceleration  to  be  gradual,  as  the  design  goal  required. 

During  tunnel  operation,  the  water  manometers  connected  to  the  instrumented 
blades  showed  fluctuations  in  the  pressures.  The  suction  side  of  the  blade  had  the  largest 
fluctuations.  The  occurrence  of  fluctuations  raised  the  question  of  flow  separation  on  the 
suction  side.  Additional  measurements  were  needed  to  determine  whether  the  flow  had 
separated. 


30 


c. 


FIVE-HOLE  PRESSURE  PROBE  MEASUREMENTS 


1.  Manual  Loss  and  AVDR  Calculations 

The  manual  system  was  used  first  to  make  a  downstream  survey  of  1 54  mm  width, 
with  an  interval  of  2.54  mm.  A  similar  upstream  survey  was  performed.  The  loss  and 
AVDR  were  calculated  using  the  equations  discussed  in  the  previous  chapter  and  found  in 
Appendix  A.  The  AVDR  was  1.023  and  the  loss  coefficient  was  0.029  using  the  Prandtl 
probe  as  a  reference.  The  calculated  loss  coefficient  was  similar  to  Gelder’s  experimental 
result  of  0.030  [Ref  1], 

Figure  14  shows  the  results  of  the  survey.  The  static  pressure  coefficient  showed  a 
rise  downstream  of  the  test  section  of  approximately  1%.  The  blade  wake  was  well-mixed 
out  at  the  downstream  station.  The  upstream  flow  angle  was  approximately  35.5  degrees 
and  the  downstream  flow  angle  was  approximately  2.5  degrees.  The  nondimensional 
velocity  showed  a  decrease  from  approximately  0.10  to  0.08.  These  experimental  results 
are  slightly  different  from  the  LDV  results  since  the  measuring  station  was  located 
approximately  2  blade  chords  upstream  and  downstream. 

2.  HP  Automated  Data  Acquisition  System 

The  HP  automated  data  acquisition  system  results  which  were  calculated  by 
Armstrong’s  “LOSS”  program  were  inconclusive.  The  calculated  loss  values  ranged  from 
-0.013  to  0.011.  The  AVDR  values  ranged  from  1.04  to  1.05. 

During  the  surveys,  the  digital  voltmeter  showed  drift  in  both  the  zero  and  span 
calibration  readings.  Additionally,  data  sampling  was  observed  to  occur  prior  to  the 
pressure  reading  from  the  Scanivalve  transducer  stabilizing  after  a  step  to  the  following 
channel.  The  conclusion,  following  these  observations  and  from  the  resulting  negative 
losses,  was  that  the  HP  system  results  were  invalid. 
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D. 


LDV  RESULTS 


1.  Inlet  Surveys 

LDV  measurements  upstream  of  the  test  section  were  performed  at  Stations  1,  2 
and  3.  Stations  1  and  3  will  be  examined  here  in  order  to  characterize  test  section  inlet 
flow  in  both  the  near  and  far  field.  Station  1  was  located  upstream  of  the  test  section  at 
30%  axial  chord  (0.30cac).  Station  1  was  surveyed  over  167%  of  passage  width,  or  254 
mm.  Three  thousand  data  points  were  taken  at  each  position  of  the  survey,  with  a  total  of 
41  positions  spaced  6.35  mm  apart.  Results  at  Station  1  in  the  form  of  velocity  ratios 
referenced  to  the  inlet  velocity  condition,  Vref,  turbulence  percentage  referenced  to  Vref, 
and  the  Reynolds  stress  correlation  coefficient,  Cuv,  are  plotted  in  Figure  15. 

Station  1  total  velocity  ratio,  WA^ref,  was  nearly  uniform  across  the  passage  span. 
Both  the  axial  velocity  ratio,  UA^ref,  and  the  tangential  velocity  ratio,  V/Vref,  showed  a 
slight  variation  across  the  passage.  The  potential  influence  of  the  blades  was  felt  as  far 
upstream  as  30%  axial  chord,  which  resulted  in  the  depressions  in  velocity  spaced  one 
blade  passage  width  apart.  The  axial  turbulence,  T„,  was  measured  to  be 
uncharacteristically  high  for  the  position  in  the  flow,  while  the  tangential  turbulence.  Tv, 
was  consistent  at  2%.  The  most  likely  reason  for  the  high  turbulence  of  Tu  was  due  to  the 
laser  optics  misalignment,  which  was  resolved  during  later  surveys.  The  Reynolds-stress 
correlation  coefficient  was  below  0.1,  showing  the  flow  to  be  random  or  uncorrelated. 

The  survey  was  repeated  at  Station  1  over  a  165.1  mm  width  at  27  positions  with 
6.35  mm  interval  spacing,  or  108%  of  passage  width  due  to  the  turbulence  discrepancy 
noted  above,  with  the  results  shown  in  Figure  16.  The  mean  velocities  showed  as  much 
variation  as  was  recorded  previously.  Both  turbulence  percentages  were  now  of 
equivalent  magnitude,  indicating  that  the  optics  problem  from  the  first  survey  was 
resolved.  The  correlation  coefficient  was  again  less  than  0. 1 . 

Station  3  was  located  upstream  5%  of  an  axial  chord  (0.05cac)  from  the  leading- 
edge.  Two  surveys  were  performed.  For  the  first  survey  with  the  laser  horizontal,  the 
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survey  width  was  only  91%  of  blade  passage  width  or  142  mm,  and  72  data  points  spaced 
2  mm  apart  were  taken.  The  survey  did  not  cover  the  full  passage  due  to  blade 
interference  with  the  two  laser  beams.  The  survey  was  repeated  at  Station  3  over  a  width 
of  203.2  mm,  with  41  positions  at  a  spacing  of  5.08  mm,  with  the  laser  pitched  upward  5 
degrees.  Figure  1 7  shows  the  overlay  of  the  data  of  the  two  surveys,  with  and  without 
laser  pitch.  The  figure  shows  the  repeatability  of  the  velocity,  turbulence  and  correlation- 
coefficient  data  for  the  two  surveys  was  excellent.  The  proximity  of  the  leading  edge  of 
the  blades  is  evident  in  all  three  plots,  shown  by  the  velocity  gradients  near  y/S  =  0  and  1 . 
The  total  velocity  decreased  to  a  minimum  along  the  stream  line  leading  to  the  stagnation 
point  on  the  leading-edge  of  the  blade,  and  increased  dramatically  on  either  side  as  the 
flow  proceeded  around  the  leading-edge.  The  maximum  total  velocity  was  to  the  right  of 
the  stagnation  point,  indicating  the  suction  side  of  the  blade.  Periodicity  was  also  shown 
by  the  repetition  of  the  velocity  distributions.  Turbulence  remained  constant  at 
approximately  2%,  with  slight  increases  near  the  stagnation  streamlines.  The  correlation 
coefficient  showed  a  peak  near  the  stagnation  points  of  the  leading-edge  on  the  pressure 
side  of  the  blades,  and  remained  less  than  0.2  throughout  the  survey,  indicating  a  slight 
reorientation  of  the  turbulence  around  the  leading-edge. 

2.  Passage  Surveys 

Passage  surveys  were  conducted  at  Stations  4  through  10  in  the  passage  between 
blades  3  and  4  with  neither  laser  yaw  nor  laser  pitch  used  during  these  surveys.  Flow  at 
Stations  5,  7,  8,  9  and  10  will  be  discussed  in  this  section,  and  all  reduced  passage  data 
can  be  found  in  Appendix  C. 

Figure  18  shows  the  Station  5  results.  The  velocity  ratios  were  smooth  and 
showed  a  decrease  across  the  passage  from  the  suction  side  to  the  pressure  side  of  the 
passage.  The  turbulence  percentages  were  nearly  equal  and  maintained  a  value  of 
approximately  2%.  The  correlation  coefficient  ranged  from  approximately  -0. 1  to  0.2. 
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The  Station  7  survey  results  are  presented  in  Figure  19.  The  velocity  data  in  the 
figure  showed  a  decrease  for  two  points  to  the  left  of  the  maximum  in  both  total  and  axial 
velocity  on  the  suction  side  of  the  blade.  The  velocities  then  decreased  gradually  as  the 
distance  away  from  the  suction  surface  increased.  Turbulence  peaked  at  12%  in  the  axial 
direction  and  decreased  to  2%.  Turbulence  remained  at  2%  in  the  tangential  direction. 
The  correlation  coefficient  ranged  from  0  to  0.2.  These  results  indicated  that  the  passage 
survey  for  Station  7  reached  slightly  into  the  boundary  layer,  which  suggested  that  the 
boundary  layer  was  growing  in  thickness,  compared  with  Station  5. 

Station  8  results  are  presented  in  Figure  20.  The  velocity  ratios  were  smooth,  with 
a  slight  decrease  in  magnitude  as  y/S  increased.  Turbulence  was  approximately  2%, 
similar  to  previously  discussed  stations.  The  correlation  coefficient  ranged  from  0.1  to 
0.3. 

Figure  21  shows  the  results  from  the  survey  at  Station  9.  The  velocity  ratios 
remained  smooth,  turbulence  remained  at  2%,  and  the  correlation  coefficient  again  ranged 
from  approximately  0.1  to  0.3.  Figure  22  shows  Station  10  results.  The  velocity  ratios 
were  similar  to  those  at  Station  9,  the  turbulence  remained  approximately  2%,  and  the 
correlation  coefficient  ranged  from  0.1  to  0.3,  as  before. 

3.  Wake  Surveys 

Wake  surveys  were  performed  at  Stations  11,  12  and  13.  The  results  from 
Stations  11  and  13,  the  near  and  far  wake,  will  be  discussed.  A  series  of  three  surveys 
were  performed  at  each  station,  ranging  in  extent  from  a  full  passage  width  and  decreasing 
to  the  width  of  just  the  wake. 

Figure  23  shows  Station  1 1  results  for  nearly  the  full  passage  width.  The  laser 
remained  horizontal  for  this  survey,  so  blade  interference  with  the  laser  beams  prevented  a 
full  passage-width  survey.  The  velocity  ratios  are  shown  to  be  nearly  uniform. 
Turbulence  remained  approximately  2%  until  the  wake  was  reached,  at  which  point  the 
turbulence  began  to  increase.  The  correlation  coefficient  remained  below  0.3. 
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The  laser  was  pitched  5  degrees  down  for  the  next  two  wake  surveys  at  Station 
1 1 .  Figures  24  and  25  show  the  results.  The  total  velocity  ratio  decreased  and  then 
increased  again  as  the  wake  was  traversed  during  the  survey.  The  axial  velocity  ratio  was 
similar  to  the  total  velocity  ratio.  The  tangential  velocity  ratio  showed  an  initial  increase 
and  then  decreased  below  zero  as  the  wake  was  surveyed.  The  wake  turbulence  showed 
two  distinct  peaks  in  the  axial  direction  as  the  wake  was  traversed,  while  the  tangential 
direction  showed  a  single  peak.  The  maximum  axial  turbulence  was  18%,  and  the 
maximum  tangential  turbulence  was  16%.  The  correlation  coefficient  started  at  a 
magnitude  of  0.2  and  became  -0.2  as  the  wake  was  traversed,  and  returned  to  a  value  of 
approximately  0.  Figure  26  shows  the  results  of  the  two  Station  1 1  wake  surveys  plotted 
together.  The  data  of  the  two  surveys  correlated  very  well,  as  shown  by  the  close  overlay. 

Station  13  was  surveyed  over  three  different  ranges  from  an  entire  blade  passage 
width  to  just  the  width  of  the  wake.  Figure  27  shows  the  first  survey  over  more  than  the 
passage  width.  The  velocity  ratio  showed  a  decrease  through  the  wake,  giving  an 
indication  of  the  necessary  width  for  the  following  survey.  The  turbulence  was 
approximately  2%  over  the  passage,  until  the  wake  was  encountered.  In  the  wake, 
turbulence  increased  to  approximately  14%.  The  correlation  coefficient  decreased  from 
0.2  to  -0.2  through  the  wake,  and  remained  approximately  0.1  outside  of  the  wake. 

Figures  28  and  29  show  more  detailed  wake  surveys  at  Station  13.  The  velocity 
ratios  again  showed  a  decrease  as  the  wake  was  traversed,  and  then  an  increase.  The  axial 
turbulence  in  the  wake  showed  two  peaks,  similar  to  the  Station  1 1  results,  and  had  a 
maximum  magnitude  of  14%.  The  tangential  turbulence  showed  a  single  maximum  peak, 
similar  to  the  Station  11  results,  with  a  maximum  magnitude  of  12%.  The  correlation 
coefficient  went  from  0.2  to  -0.2  in  the  wake,  while  remaining  at  approximately  0.1 
outside  of  the  wake.  Figure  30  shows  the  two  Station  13  wake  surveys  plotted  together. 
The  data  correlated  well  and  only  the  axial  turbulence,  T„,  showed  a  slight  difference  in 
value  on  the  right  side  of  the  wake. 
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The  axial  turbulence  results  in  Figures  26  and  30  are  of  interest.  While  both 
stations  show  two  peaks  in  Tu,  Station  1 1  showed  the  maximum  peak  to  be  located  on  the 
left  side  of  the  wake,  while  Station  13  showed  the  maximum  magnitude  to  be  located  on 
the  right  hand  side.  While  asymmetry  in  the  wake  turbulence  was  measured  in  the  wake  of 
the  Stator  67A  cascade  [Ref  5],  a  reversal  of  the  asymmetry  from  one  side  to  another  was 
not  previously  observed.  ' 

4.  Boundary  Layer  Surveys 

Boundary  layer  surveys  were  performed  at  four  stations,  normal  to  the  surface  of 
the  blade  at  the  specified  station.  Station  5,  7  and  9  were  surveyed  on  the  suction  side, 
while  Station  8  was  surveyed  on  the  pressure  side.  A  combination  of  laser  pitch  and  yaw 
was  used  in  order  to  position  the  LDV  probe  volume  as  close  to  the  blade  as  possible. 

A  total  of  three  surveys  were  performed  at  Station  5  with  three  different 
combinations  of  laser  yaw  and  pitch.  The  combination  which  worked  best  was  a  yaw  of  5 
degrees  and  a  pitch  of  0  degrees,  due  to  beam  reflection  from  the  cascade  wall  and  off  the 
brass  shims  on  the  far  endwall.  The  results  of  the  boundary-layer  survey  at  Station  5  are 
plotted  in  terms  of  the  ratio  of  distance  from  the  wall  and  blade  chord,  and  are  shown  in 
Figure  31.  The  velocity  ratios  showed  only  a  slight  decrease  on  the  left  side,  at  the 
minimum  y/S  value,  which  corresponded  to  the  edge  of  the  boundary  layer.  The 
turbulence  was  a  maximum  for  approximately  4  survey  points  at  the  boundary  layer  edge, 
which  decreased  to  approximately  2%  away  from  the  boundary  layer.  The  correlation 
coefficient  was  a  maximum  near  the  blade  surface,  and  decreased  to  0. 1 .  No  data  could  be 
taken  close  to  the  blade. 

The  results  of  the  boundary-layer  survey  taken  at  Station  7  are  shown  in  Figure  32. 
The  velocity  ratios  showed  a  low  flat  area  near  the  blade,  which  then  rose  to  near 
freestream  values.  The  axial  turbulence  started  near  12%,  rose  to  a  peak  of  40%  and  then 
began  to  decrease  near  the  freestream.  Tangential  turbulence  remained  approximately  2%. 
The  correlation  coefficient  fluctuated  from  -0. 1  to  0. 1  for  the  survey.  The  high  level  of 
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axial  turbulence  suggested  that  the  flow  had  separated  in  the  vicinity  of  Station  7.  The 
boundary  layer  data,  along  with  the  fluctuating  pressure  noted  during  Cp  measurements, 
and  the  passage  survey  results,  collectively  indicated  flow  separation. 

The  results  from  the  Station  8  pressure-side  boundary-layer  survey  are  shown  in 
Figure  33.  The  abscissa  was  plotted  in  reverse  order,  since  the  pressure  side  of  blade  4 
was  on  the  right-hand  side  of  the  survey  passage.  The  velocity  ratios  showed  a  decrease 
in  the  total  velocity  for  the  last  seven  points  of  the  survey,  on  the  right  side.  The 
turbulence  increased  for  the  seven  positions  near  the  blade,  to  a  maximum  of  6%.  The 
correlation  coefficient  showed  an  increase  from  approximately  0.3  to  0.4  in  the  boundary 
layer. 

Station  9  boundary  layer  results  are  plotted  in  Figure  34.  This  survey  was 
performed  on  the  suction  side  of  the  blade.  The  velocity  ratio  showed  a  decrease  for  the 
ten  points  near  the  blade  surface,  and  the  turbulence  increased  in  both  the  axial  and 
tangential  direction  as  the  blade  surface  was  approached.  The  boundary  layer  profile 
suggested  that  the  flow  had  reattached  itself  to  the  suction  surface  prior  to  Station  9.  The 
correlation  coefficient  decreased  below  0  in  the  boundary  layer,  and  then  increased  to  0. 1 
in  the  passage. 


5.  Combination  of  Boundary  Layer  and  Passage  Surveys 

Two  Station  7  survey  results  are  presented  together  in  Figure  35.  The  passage 
survey  at  Station  7  indicated  the  boundary  layer  was  thicker  than  at  other  stations, 
allowing  a  few  data  points  to  be  taken  inside  the  boundary  layer  without  laser  pitch  or 
yaw.  The  boundary  layer  survey  at  Station  7,  performed  along  a  path  normal  to  the  blade 
surface  at  Station  7,  showed  that  a  large  region  of  separation  existed.  In  order  to  study 
the  flow  close  to  the  blade  surface  with  the  same  reference  as  the  passage  surveys,  the 
laser  was  yawed  to  perform  another  passage  survey  at  Station  7.  Data  from  this  survey 
with  the  laser  yawed  was  plotted  with  the  Station  7  passage  survey  data  [from  Figure  19], 
in  Figure  35.  The  velocity  data  in  the  figure  show  a  flat  region  for  both  total  and  axial 
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velocity  in  the  boundary  layer,  and  then  a  rapid  rise  to  freestream  conditions.  Turbulence 
peaked  at  24%  in  the  axial  direction,  and  remained  at  2%  in  the  tangential  direction.  The 
correlation  coefficient  was  erratic  in  the  boundary  layer,  ranging  from  -0. 1  to  0. 1 .  The 
correlation  coefficient  then  settled  to  a  positive  value  of  0.2  outside  the  boundary  layer. 

Hobson  and  Shreeve  [Ref  9]  located  separation  near  the  leading  edge  of  Stator 
67A  blades  in  the  NPS  cascade.  The  velocity  profile  at  Station  3  from  Hobson  and 
Shreeve  [Ref  9]  had  a  similar  shape  to  the  Stator  67B  Station  7  velocity  profile  [Figure 
35].  The  turbulence  results  from  their  measurements  and  those  at  Station  7  are  quite 
similar.  The  results  of  the  combined  data  in  Figure  35,  Cp  observations,  and  the  boundary 
layer  data  all  indicate  that  flow  separation  was  present  on  the  suction  side  at  Station  7. 

The  passage  survey  data  and  the  boundary  layer  data  at  Stations  8  and  9  are 
plotted  in  Figures  36  and  37,  respectively.  The  results  were  of  interest  since  the  passage 
survey  was  nearly  normal  to  the  blade  at  those  stations.  The  boundary  layer  extent  was 
evident  in  these  figures.  The  velocity  ratios  decreased  in  the  boundary  layer  and  the 
turbulence  increased  significantly. 


E.  CFD  RESULTS 

Rotor  Viscous  Code  Quasi-3-D  was  nin  using  three  turbulence  models.  The 
pressure  ratio,  prat,  was  iterated  upon  until  the  inlet-flow  angle  converged  to  within  0.02 
degrees.  The  CFD  predictions  are  plotted  in  Figure  38  together  with  the  experimental 
results  for  the  pressure  distribution.  Table  3  gives  the  values  of  the  pressure  ratio,  prat, 
which  yielded  design  inlet-flow  angle  when  the  CFD  code  was  used  with  the  three 
different  turbulence  models. 

All  computed  results  showed  a  sharp  acceleration  on  the  suction  side  initially,  and 
then  a  gradual  deceleration  for  most  of  the  rest  of  the  suction  side.  Near  ^/c  =  0.70  on  the 
suction  side,  the  code  predicted  a  slight  re-acceleration.  These  CFD  results  for  the  suction 
side  did  not  match  the  experimental  data  which  showed  a  more  gradual  acceleration. 
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peaking  at  approximately  ^/c  =  0.40,  and  then  a  near-linear  deceleration.  The  predicted 
results  on  the  pressure  side  of  the  blade  compared  quite  well  with  the  experimental  data. 

The  Wilcox  k-omega  turbulence  model  had  the  capability  to  modify  the  freestream 
turbulence,  which  was  computationally  changed  from  0.02  to  0.04.  The  result  from  an 
increase  in  turbulence  from  0.02  to  0.04  showed  no  significant  change  in  the  predicted  Cp 
results. 
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Figure  14.  Manual  Loss  Survey  Results. 
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Figure  16.  Repeat  Station  1  Survey  Results  for  108%  Passage  Width. 
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Figure  17.  Station  3  Survey  Results  With  and  Without  Laser  Pitch. 
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Figure  20.  Station  8  Passage  Survey  Results. 
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Figure  22.  Station  10  Passage  Survey  Results. 
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Figure  23.  Station  11  Survey  Results. 
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Figure  24.  Station  1 1  Wake  Surv^  Number  1  Results. 
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Figure  26.  Station  1 1  Wake  Surveys  1  and  2  Results. 
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Figure  29.  Station  13  Wake  Survey  Number  2  Results. 
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Figure  30.  Station  13  Wake  Surveys  1  and  2  Results 


Figure  3 1 .  Station  5  Boundary  Layer  Survey  Results. 
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Figure  33.  Station  8  Boundary  Layer  Survey  Results. 
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Vefocity  Ratio 


Figure  34.  Station  9  Boundary  Layer  Survey  Results. 
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station  7  Velocity  Ratios 


Figure  35.  Station  7  Boundary  Layer  and  Passage  Survey  Results. 
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station  8  Velocity  Ratios 


Figure  36.  Station  8  Boundary  Layer  and  Passage  Survey  Results. 


Correlation  Coefficient  Turbulence  (%)  Velocity  Ratio 


Model 

Cebeci-Smith 

Baldwin-Lomax 

Wilcox  k-omega 

Pressure  Ratio  (prat) 

0.9765 

0.97635 

0.97625 

Table  3.  Pressure  Ratios  for  Turbulence  Models. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Compressor  Stator  67B  cascade  blading  was  successfully  installed  and 
experimentally  tested  at  design  conditions  in  the  NPS  Low-Speed  Cascade  Wind  Tunnel. 
The  experiments  were  conducted  at  a  design  inlet-flow  angle  of  36.3  degrees,  a  Mach 
number  of  0.22  and  a  Reynolds  number  based  on  chord  of  640,000. 

Experimental  blade  surface  pressure  measurements  were  obtained  using  water 
manometers  and  an  HP  automated  data  acquisition  system.  The  resulting  Cp  plots 
showed  that  the  flow  on  the  suction  side  rapidly  accelerated  from  the  leading  edge,  more 
gradually  accelerated  until  midspan,  and  then  gradually  decelerated  to  the  trailing  edge. 
Fluctuations  in  pressure  on  the  suction  side  of  the  blades  observed  during  the  experiments 
suggested  unsteady  flow  or  separation  had  occurred  on  the  blade,  and  this  was  later 
confirmed  with  LDV  passage  and  boundary  layer  surveys. 

Experimental  measurements  yielded  a  loss  of  0.029.  These  results  compared  well 
to  NASA  LeRC  loss  values  of  0.030  [Ref  1].  However,  the  NASA  results  were  from 
tests  of  the  entire  stage,  so  the  comparison  is  not  totally  valid. 

The  LDV  results  characterized  the  flow  in  the  inlet,  in  the  blade  passage,  in  the 
boundary  layers  and  at  the  exit  of  the  test  section  at  design  conditions.  Inlet  surveys 
measured  the  upstream  influence  of  the  blades.  Passage  surveys  showed  a  smooth 
acceleration  up  to  a  separation  point  in  the  vicinity  of  Station  7  on  the  suction  side  of  the 
blade,  and  then  a  smooth  deceleration  and  reattachment  of  the  flow.  With  the  presence  of 
separation,  although  the  blade  produced  a  useful  force,  the  controlled  diffusion  design 
goal  of  preventing  separation  at  design  incidence  was  not  achieved.  Boundary-layer 
surveys  were  successfully  performed  at  three  suction-side  stations,  and  one  pressure-side 
station.  The  boundary-layer  surveys  defined  the  growth  of  the  boundary  layer  throughout 
the  passage.  Wake  surveys  showed  the  influence  of  mixing  in  the  wake. 
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CFD  predictions  of  the  blade  surface  pressure  coefficients  using  RVCQ3D  Version 
300  did  not  correlate  well  with  the  experimental  results.  Three  different  turbulence 
models  were  used,  with  similar  results. 


B.  RECOMMENDATIONS 

More  extensive  boundary-layer  surveys  and  analysis  should  be  performed  at  design 
incidence  to  determine  more  exactly  the  regions  of  separation  and  reverse  flow.  Laser 
rotation,  in  addition  to  laser  pitch  and  yaw,  could  allow  the  LDV  probe  volume  to  be 
positioned  closer  to  the  blade  surface  at  the  measurement  stations.  Flow  visualization 
using  a  laser  sheet  should  be  performed  to  help  determine  the  separation  location,  and  to 
determine  the  extent  of  the  separation  region. 

The  five-hole  probe  data  acquisition  system  should  be  improved  by  updating  the 
survey  equipment,  and  controlling  software.  A  computer-controlled  automatic-traverse 
mechanism  for  the  five-hole  probe,  and  an  automated  yaw  system  to  align  the  five-hole 
probe  with  the  flow,  would  permit  surveys  to  be  accomplished  more  rapidly.  An  update  of 
the  HP  BASIC  software  to  include  self-adjusting  time  pauses  to  allow  settling  of  pressure 
data  prior  to  recording,  point-by-point  data  recording  in  case  of  system  failure,  and  an 
option  to  use  Prandtl  probe  pressures  for  reference  should  be  implemented  immediately. 
The  option  of  using  software  based  on  “Lab VIEW”  to  control  the  system  should  be 
investigated.  Finally,  the  five-hole  probe  calibration  coefficients  should  be  verified  before 
further  probe  surveys  are  conducted. 

Additional  analysis  with  RCVQ3D  should  be  performed  to  try  to  reproduce  the 
experimental  results.  Different  grid  sizes  and  number  of  iterations  should  be  investigated. 
Also,  the  incidence  angle  should  be  varied  in  an  attempt  to  match  the  inlet  flow  conditions. 
Loss  calculations  should  be  included  in  the  CFD  analysis.  Experimentation  with  a  full 
three  dimensional  code  might  provide  better  results. 
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When  the  flow  at  the  design  incidence  angle  is  initially  characterized,  positive  off- 
design  incidence  angles  should  be  set  to  determine  immediately  the  available  range  of 
incidence  and  the  on-set  of  stall.  This  will  allow  the  future  test  plan  to  be  determined  most 
optimally. 
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APPENDIX  A.  FIVE-HOLE  PROBE  EQUATIONS 


Dimensionless  velocity  coefficient:  P  =  ^ 

Pi 

Average  yaw  pressure:  P23  =  -- 

2 

Pitch  angle  coefficient:  — — 

P1-P23 


Dimensionless  velocity  coefficients  [Ref  8]  Cy: 


Cii 

Ci2 

Ci3 

Ci4 

Ci5 

C.j 

0.015926 

4.932133 

-153.66876 

3137.9614 

-24299.005 

C2j 

-0.003563 

0.699505 

-62.977261 

2068.3721 

-20872.148 

C3j 

0.080098 

-24.844173 

1980.4954 

1980.4954 

541835.5 

Dimensionless  velocity  polynomial: 

x  =  {c„+c„.p+c„  (i“+c„-p’+c„  r} 

+  ^C21  +  C22  ■  p  +  C23  'P  +f'24'P  ■*'^25'P  I’T 5-hole 
+  {C3,  +C32  P  +  Cjj  "P  +Cj4-P  +C3J-P  I'ys^hole 


Reference  dimensionless  velocity: 


Reference  flow  function: 


Axial-Velocity-Density  Ratio: 


K  = 


WL 


I-X^ 


l-X: 


ref 


Y-1 


CO 


5- hole) 


AVDR  = 


fK.dx 
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Dimensionless  total  pressure: 


Dimensionless  static  pressure: 


Loss  coefficient: 


P. 

P. 


c  — 

p 


o  = 


f  - ^ —  f 

Jo  AVDR  ^ 


jl  Cp,^K„dx  -  JCp,„3K„,dx 
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APPENDIX  C.  LDV  SUMMARY  AND  REDUCED  DATA 
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APPENDIX  D.  REFERENCE  VELOCITY  CODE 


FORTRAN  CODE  “CALIBl” 


C 


C 

PROGRAM  TO  COMPUTE  THE  CALIBRATION  CURVE  FOR  THE  NPS  LOW  SPEED  C 
CASCADE  WIND  TUNNEL. 

C 

C 

- C 


A  STRAIGHT  LINE  IS  FITTED  THROUGH  THE  REFERENCE  CONDITIONS  OF 
THE  TUNNEL  AT  DIFFERENT  SPEEDS. 

THE  REFERENCE  VELOCITY  IS  THEN  OBTAINED,  BY  NEWTON'S  METHOD, 
DEPENDING  ON  THE  TUNNEL  PLENUM  PRESSURE  AND  TEMPERATURE 


- c 

PROGRAM  CALIBRATE 
IMPLICIT  REAL*8(A-H,0-Z) 

PARAMETER  (NP=6) 

DIMENSION  VA(NP),VT(NP),PA(NP),PP(NP),TP(NP) 

DIMENSION  VTOT(NP),X(NP),ANUX(NP),PR(NP) 

DIMENSION  PAR(100),PPR{100),TPR(100),PRR(100),ANUXR(100) 

DIMENSION  VREF(IOO) 

CHARACTER*  14  NAME(IOO) 

OPEN(UNIT=10,FILE='CALIB.DAT',STATUS='UNKNOWN’) 

OPEN(UNIT=ll,FILE=’REFER.DAr,STATUS='UNKNOWN') 

OPEN(UNIT=12,FILE='CALIB.OUT',STATUS=’UNKNOWN') 

PRINT  BANNER 

WRITE(*,500) 

WRITE(12,500) 

DOFORMAT(1X,78('C'),/1X,'C,76X,’C/, 

#1X,'C',20X, 'OUTPUT  FROM  PROGRAM  CALIBRATE',27X,'C'/1X,'C',76X,'C'/ 
#1X,78('C')// 

#5X,'LEAST  SQUARES  STRAIGHT  LINE  CURVE  FIT  IS  USED'/ 

#5X,TO  DETERMINE  TUNNEL  CHARACTERISTICS  AT  DIFFERENT  SPEEDS'// 

#5X, 'NEWTON  S  METHOD  IS  USED  TO  DETERMINE  THE  REFERENCE  VELOCITY'/ 
#5X,'FROM  THE  RECORDED  AMBIENT  PRESSURE  AND  TUNNEL  PLENUM'/ 
#5X,'PRESSURE  AND  TEMPERATURE'/) 

INITIALIZE  AIR,  MERCURY  AND  WATER  PROPERTIES 


105 


n  n  n  n  n  n  n 


c 

RHOW=1000.D0 

RHOHG=13000.D0 

RHOA=1.2256DO 

CPA=1005.D0 

GAMMA=1.4D0 

GM1=GAMMA-1.D0 

C 

C  BEGIN  DETERMINING  TUNNEL  CHARACTERISTICS 
C 

WRITE(*50I) 

WRITE(I2,501) 

501  F0RMAT(//1X,’BEGIN  DETERMINING  TUNNEL  CHARACTERISTICS',/ 
#1X,'FR0M  THE  FOLLOWING  MEASURED  VALUES’,// 

#3X,'AXIAL  VEL.  TANGENTIAL  VEL.  AMBIENT  PRESS.' 

#3X,’PLENUM  PRESS.  PLENUM  TEMP. '/ 

#3X,'M  PER  SEC.  M  PER  SEC.  INCHES  MERCURY' 

#3X,'INCHES  WATER  DEG.  C.’/) 

C 

C  READ  IN  DATA  POINTS 
C 

C  VA  =  AXIAL  VELOCITY  -  MEASURED  BY  LDV  (M/S) 

C  VT  =  TANGENTIAL  VELOCITY  -  MEASURED  BY  LDV  (M/S) 

C  PA  =  AMBIENT  PRESSURE  (INCHES  MERCURY) 

C  PP  =  PLENUM  PRESSURE  (INCHES  WATER) 

C  TP  =  PLENUM  STAGNATION  TEMPERATURE  (DEG.  C.) 

C 

DO  1  I=1,NP 

READ(10,100)VA(I),VT(I),PA(I),PP(I),TP(I) 

WRITE(*,  1 1 0)  VA(I),  VT(I),P  A(I),PP(I),TP(I) 

WRITE(I2, 1 10)VA(I),VT(I),PA(I),PP(I),TP(I) 

1  CONTINUE 

100  F0RMAT(1X,5F8.4) 

1 10  FORMAT(3X,F8.4,7X,F8.4,8X,F8.4,7X,F8.4,7X,F8.4) 

C 

C  CALC  VTOT,  X,  ANUX  AND  PR 
C 

WRITE(*,510) 

WRITE(12,510) 

510  FORMAT(/5X,'CALCULATED  VALUES  FOR  THE  TUNNEL  CONnGURATION',// 
#3X,TOTAL  VELOCITY’,6X,’MACH  NUMBER',9X,'MACH  NUMBER  FUNCT.’,3X, 
^’PRESSURE  RATIO'//) 

C 

DO  2  I=1,NP 

VTOT(I)=DSQRT(VA(I)**2+VT(I)**2) 

X(I)=VTOT(I)/DSQRT(2.DO*CPA*(TP(I)+273.16DO)) 

ANUX(I)=(GAMMA/GM1)*(X(1)**2)*(LD0-(X(I)**2))**(1.D0/(GM1)) 

PR(I)=  I  .DO-(RHOHG*PA(I))/(RHOW*PP(I)+RHOHG*PA(I)) 

WRITE(*,  1 0 1 )  VTOT(I),X(I),  ANUX(I),PR(I) 

WRITE!  12,101 )  VTOT(I),X(I),  ANUX(I),PR(I) 

2  CONTINUE 
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101  FORMAT(4D20.11) 

C 

C  CALL  THE  LEAST  SQUARES  SUBROUTINE  TO  FIT  A  STRAIGHT  LINE  THROUGH  THE 
C  DATA;  X-AXIS  =  PR,  Y-AXIS  =  ANUX  (MACH  NO,  PARAMETER) 

C 

WRITE(*,520) 

WRITE(12,520) 

520  FORMAT(/10X, 'CALLING  LEAST  SQUARES  SUBROUTINE', /IX, 

#'TO  DETERMINE  THE  PRESSURE  RATIO  AS  A  FUNCTION  OF  MACH  NO.  PARAM' 
#//lX, 'PRESSURE  RATIO  =  A1  ♦  ANUX  +  AO',//) 

C 

CALL  LEASTSQUARE(NP,ANUX,PR,A0,A1) 

C 

WRITE(*,530)A1,AO 

WRITE(12,530)A1,A0 

530  FORMAT(/10X,'A1  =  ',D20. 11,’  AO  =  ',020. 1  If) 

C 

C  READ  IN  REFERENCE  CONDITIONS 
C 

WRITE(*,199) 

WRITE(12,199) 

199  FORMAT(/lX, 'REFERENCE  CONDITIONS  FOR  EACH  RUN'/ 

#1X,'AMBIENT  PRESSURE  PLENUM  PRESSURE  PLENUM  TEMPERATURE', 
#3X,'RUNNAME'/ 

#1X,' INCHES  MERCURY  INCHES  WATER  DEGREES  CELSIUS’//) 

C 

C  THE  NUMBER  OF  EXPERIMENTS 
C 

READ(11,198)NE 
198  FORMAT(lX,I3) 

C 

C  THE  REFERENCE  CONDITIONS 
C 

DO  99  I=1,NE 

READ(11,200)PAR(I),PPR(I),TPR(I),NAME(I) 

WRITE(*,201)PAR(I),PPR(I),TPR(I),NAME(I) 

WRITE(12,201)PAR(I),PPR(I),TPR(I),NAME(I) 

99  CONTINUE 

200  F0RMAT(1X,3F8.4,A14) 

201  FORMAT(2(6X,F8.4),10X,F8.4,1 1X,A14) 

C 

C  CALCULATE  THE  PRESSURE  RATIO 
C 

DO  88  I=1,NE 

PRR(I)=  1  .DO-(RHOHG'’'PAR(I))/(RHOW*PPR(I)+RHOHG*PAR(I)) 

C 

ANUXR(I)=(-A0+PRR(I))/A  1 
C 

WRITE(*, 550)1, PRR(I),ANUXR(I),NAME(I) 

WRITE(  1 2,550)1, PRR(I),ANUXR(I),NAME(I) 

550  FORMAT(/lX,'I  =  ’,13,/ 
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#'  PRESSURE  RATIO  =  ',F10.5,’  MACH  NUMBER  PARAMETER  =  •,E12.4,/ 

#'  RUN  NAME  =  ’,A14/) 

C 

C  NEWTON  METHOD  TO  DETERMINE  THE  ROOTS  OF  THE  EQUATION  FOR  ANUX 
C 

C  GUESS  INITIAL  VALUE  OF  VREF  =  VTOT(5) 

C 

WRITE(*,560) 

WRITE(12,560) 

560  FORMAT(//IOX, 'BEGIN  NEWTON  ITERATION'/) 

C 

KOUNT=l 

VREF(I)=VTOT(5) 

XR=VREF(I)/DSQRT(2.D0*CPA*(TPR(I)+273.16D0)) 

1000F=(GAMMA/GMI)*(XR**2)*(1.D0-(XR**2))**(1.D0/(GM1))-ANUXR(I) 

DFDX=(GAMMA/GM1)*2.D0'»XR*((LD0)-(XR'»'»2)'»(1.D0/(GMI)))* 

#  ( I  .D0-(XR'^*2)*(I  .DO/GM I  )*(I  .D0/(1  .D0-(XR'»'»2)))) 

TERM=F/DFDX 

C 

WRITE('»,570)KOUNT,XR,TERM 

WRITE(12,570)KOUNT,XR,TERM 

570  FORMAT(lX,'ITERATION  NUMBER  '.12,'  MACH  NO.  PARAM.  =  ',F8.6, 

#'  ERROR  TERM  =  ',D  1 2.4) 

C 

XR=XR-(F/DFDX) 

KOUNT=KOUNT+I 
IF(DABS(TERM).LT.1.D-10)GO  TO  999 
IF(KOUNT.LT.10.OR.DABS(TERM).GT.l.D-10)GO  TO  1000 
C 

C  CALCULATE  THE  REFERENCE  VELOCITY 
C 

999  VREF(I)=XR*DSQRT(2.D0*CPA*(TPR(I)+273. 16)) 

WRITE(*,580)VREF(I) 

WRITE(  1 2,580)  VREF(I) 

580  FORMAT(/lX,'VREF  =  ',F20. 1 1/) 

88  CONTINUE 
C 

C  PRINT  FINAL  RESULT 
C 

WRITE(  12,585) 

585  FORMAT(/lX,'EXPERIMENT  NUMBER  REFERENCE  VELOCITY  NAME'/) 

DO  777  I=1,NE 

WRITE(  12,590)1,  VREF(I),NAME(I) 

777  CONTINUE 

590FORMAT(8X,I2,16X,F8.4,10X,A14) 

STOP 

END 

C - C 

c 
c 

C  A  LEAST  SQUARES  CURVE  FIT  FOR  A  STRAIGHT  LINE  THROUGH  NOISY 
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non 


C  DATA 
C 

C  ALGORITHM  TAKEN  FROM:  NUMERICAL  METHODS,  ROBERT  W.  HORNBECK 
C  PAGES  122  -  130,  QUANTUM,  1975 

C 


SUBROUTINE  LEASTSQUARE(N,X,FX,A0,A1) 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  COEFF(2,2),RHS(2),X(N),FX(N) 

C 

C  SET  UP  COEFFICIENT  MATRIX  AND  RIGHT  HAND  SIDE. 
C 

COEFF(l.l)=N 

COEFF(1,2)=O.DO 

COEFF(2,2)=O.DO 

C 

RHS(1)=0.D0 

RHS(2)=0.D0 

C 

DO  2  I=1,N 

COEFF(  1 ,2)=COEFF(  1 ,2)+X(I) 

COEFF(2,2)=COEFF(2,2)+X(I)**2 

RHS(1)=RHS(1)+FX(I) 

RHS(2)=RHS(2)+X(I)*FX(I) 

2  CONTINUE 
COEFF(2, 1  )=COEFF(  1,2) 

C 

C  PRINT  MATRIX  EQUATION 
C 

WRITE(*,110) 

WRITE(12,110) 

110  FORMAT(/5X,'MATRIX  EQUATION',/) 

DO  3  1=1,2 

WRITE(*,102)(COEFF(LJ),J=1,2),I-1,RHS(1) 

WRITE(12,102)(COEFF(1,J),J=1,2),I-1,RHS(I) 

3  CONTINUE 

102  FORMAT(5X,2(2X,D10.2),'  A',I1,3X,D10.2) 

C 

C  GAUSS  ELIMINATION 
C 

TERM=COEFF(2,  l)/COEFF(  1, 1) 

C 

COEFF(2,2)=COEFF(2,2)-COEFF(  1 ,2)*TERM 
RHS(2)=RHS(2)-RHS(  1 )  *TERM 
C 

Al=RHS(2)/COEFF(2,2) 

A0=(RHS(  1  )-COEFF(  1 ,2)*  A 1  )/COEFF(  1,1) 

C 

RETURN 

ENDD 
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n  n  n  n  n 


FORTRAN  INPUT  FILE  “REFER.DAT’ 


28 

29.9499  12.1000  21.1111  0502sl3a 

29.9499  11.9000  21.1111  0502sl3b 

29.9499  11.9000  21.1111  0502sl3c 

29.9499  11.9000  21.1111  0502sl2a 

29.9499  11.9000  21.1111  0502sl2b 
29.8684  12.0000  20.0000  0506slla 

29.8481  12.1000  20.5556  0506sl0j 

29.8481  12.1000  20.5556  0506s9_ 

29.8481  12.1000  20.5556  0506s8_ 

29.8481  11.9000  21.6667  0506s7a_ 
30.0517  12.0000  21.1111  0509s6a_ 
30.0313  12.1000  20.5556  0509s5_ 

29.8888  12.0000  21.1111  0509s4_ 

29.8888  11.9500  20.5556  0507s3a_ 

29.8888  11.9000  20.0000  0507s2a_ 
29.9702  12.0000  21.1111  0401sla_ 

29.9295  1 1.9000  21.1111  0520s6b_ 

29.9295  11.8000  22.2222  0520sllb 

29.9295  11.9500  22.2222  0520sllc^ 
29.9091  12.2000  22.2222  0520bp8_ 

29.9499  11.9500  20.0000  0526s2b_ 

29.9499  12.0000  20.5556  0526s3b_ 

29.9295  11.9500  20.5556  0603bs5_ 
29.9091  12.0000  25.5556  0708bs9_ 

29.9295  11.9000  23.3333  0404slb_ 
29.8684  12.0000  22.2222  0901bs7_ 

29.8481  11.8500  24.4444  0901bs7b_ 

29.8481  11.8000  25.0000  0901s7b 


no 


FORTRAN  OUTPUT  FILE  “CALIB.OUT 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  OUTPUT  FROM  PROGRAM  CALIBRATE  C 

c  c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

LEAST  SQUARES  STRAIGHT  LINE  CURVE  FIT  IS  USED 
TO  DETERMINE  TUNNEL  CHARACTERISTICS  AT  DIFFERENT  SPEEDS 

NEWTON  S  METHOD  IS  USED  TO  DETERMINE  THE  REFERENCE  VELOCITY 
FROM  THE  RECORDED  AMBIENT  PRESSURE  AND  TUNNEL  PLENUM 
PRESSURE  AND  TEMPERATURE 


BEGIN  DETERMINING  TUNNEL  CHARACTERISTICS 
FROM  THE  FOLLOWING  MEASURED  VALUES 


AXIAL  VEL.  TANGENTIAL  VEL.  AMBIENT  PRESS.  PLENUM  PRESS.  PLENUM  TEMP. 
M  PER  SEC.  M  PER  SEC.  INCHES  MERCURY  INCHES  WATER  DEG.  C. 


33.4750 

24.5991 

29.9499 

3.6000 

20.5556 

42.7350 

31.4532 

29.9499 

5.8000 

21.6667 

50.4851 

37.3610 

29.9499 

8.1000 

22.7778 

56.1076 

41.5132 

29.9295 

10.1000 

23.3333 

60.7966 

45.0193 

29.9295 

11.9000 

23.3333 

65.9383 

48.5611 

29.9295 

14.0000 

23.8889 

CALCULATED  VALUES  FOR  THE  TUNNEL  CONFIGURATION 
TOTAL  VALOCITY  MACH  NUMBER  MACH  NUMBER  FUNCT.  PRESSURE  RATIO 


0.41541441306E+02 
0.53062076997E+02 
0.628059682 12E+02 
0.69795476587E+02 
0.756502739I9E-H)2 
0.81890413603E+02 


0.54065475267E-01 
0.68929132229E-01 
0.8I4334I8389E-01 
0.90411 111093E-01 
0.97995252042E-01 
0.10597930366E+00 


0.10156165233E-01 
0.16432468008E-01 
0.2282703 1824E-01 
0.2802852 1818E-01 
0.32809629573E-01 
0.38216120664E-01 


0.91615012341E-02 
0.14678018680E-01 
0.20379988663E-01 
0.25301646987E-01 
0.2967703 1599E-01 
0.34732257708E-01 


MACH  NO.  PARAM 


CALLING  LEAST  SQUARES  SUBROUTINE 
TO  DETERMINE  THE  PRESSURE  RATIO  AS  A  FUNCTION  OF 

PRESSURE  RATIO  =  A1  *  ANUX  +  AO 


MATRIX  EQUATION 
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0.60E+01  0.15E+00A0  0.13E+00 

0.15E+00  0.42E-02A1  0.38E-02 

Al=  0.91276427032E+00  AO  =  -0.26460149149E-03 


REFERENCE  CONDITIONS  FOR  EACH  RUN 

AMBIENT  PRESSURE  PLENUM  PRESSURE  PLENUM  TEMPERATURE  RUN  NAME 
INCHES  MERCURY  INCHES  WATER  DEGREES  CELSIUS 


29.9499 

12.1000 

21.1111 

0502sl3a 

29.9499 

11.9000 

21.1111 

0502sl3b 

29.9499 

11.9000 

21.1111 

0502sl3c 

29.9499 

11.9000 

21.1111 

0502sI2a 

29.9499 

11.9000 

21.1111 

0502sl2b 

29.8684 

12.0000 

20.0000 

0506s 11 

29.8481 

12.1000 

20.5556 

0506s 10 

29.8481 

12.1000 

20.5556 

0506s9 

29.8481 

12.1000 

20.5556 

0506s8 

29.8481 

11.9000 

21.6667 

0506s7 

30.0517 

12.0000 

21.1111 

0509s6a 

30.0313 

12.1000 

20.5556 

0509s5 

29.8888 

12.0000 

21.1111 

0509s4 

29.8888 

11.9500 

20.5556 

0507s3a 

29.8888 

11.9000 

20.0000 

0507s2a 

29.9702 

12.0000 

21.1111 

0401sla 

29.9295 

11.9000 

21.1111 

0520s6b 

29.9295 

11.8000 

22.2222 

0520s 11a 

29.9295 

11.9500 

22.2222 

0520sllb 

29.9091 

12.2000 

22.2222 

0520bp8 

29.9499 

11.9500 

20.0000 

0526s2b 

29.9499 

12.0000 

20.5556 

0526s3b 

29.9295 

11.9500 

20.5556 

0603bs5 

29.9091 

12.0000 

25.5556 

0708bs9 

29.9295 

11.9000 

23.3333 

0404slb 

29.8684 

12.0000 

22.2222 

0901bs7 

29.8481 

11.8500 

24.4444 

0901bs7b 

29.8481 

11.8000 

25.0000 

0901s7b 

EXPERIMENT  NUMBER 

REFERENCE  VELOCITY  NAN 

1 

75.9548 

0502sl3a_ 

2 

75.3334 

0502sl3b 

3 

75.3334 

0502sl3c 

4 

75.3334 

0502sl2a_ 

5 

75.3334 

0502sl2b 

6 

75.6033 

0506slla_ 

7 

76.0105 

0506sl0 

8 

76.0105 

0506s9 
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9 

76.0105 

0506s8 _ 

10 

75.5311 

0506s7a_ 

11 

75.5184 

0509s6a_ 

12 

75.7816 

0509s5_ 

13 

75.7209 

0509s4_ 

14 

75.4939 

0507s3a_ 

15 

75.2668 

0507s2a 

16 

75.6195 

0401sla 

17 

75.3587 

0520s6b 

18 

75.1875 

0520sllb 

19 

75.6570 

0520s 11c 

20 

76.4587 

0520bp8j 

21 

75.3466 

0526s2b 

22 

75.5733 

0526s3b 

23 

75.4433 

0603bs5 

24 

76.2651 

0708bs9_ 

25 

75.6427 

0404s lb 

26 

75.8893 

0901bs7 

27 

75.7288 

0901bs7b 

28 

75.6417 

0901s7b  ' 
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APPENDIX  E.  GRAPE  AND  RVCQ3D  INPUT  AND  PCP  CODE 


“GRAPE”  INPUT 

&GRID1 

JMAX=340,KMAX=49,NTETYP=3,NAIRF=5,JAIRF=343,NIBDST=7, 

DSI=0.0002,  JTEBOT=80,  JTETOP=26 1  ,XLE=0.0,NOBSHP=7,XTE=  1 .0, 

ALAMF=0.0,ALAMR=0.0,XLEFT=-.5,XRIGHT=3., 

RCORN=.10,NOUT=4,NORDA=5,3,MAXITA=0,300, 

&END 

&GRID2 

AAAI=3,BBBI=3,DSOBI=0.01,ROTANG=-16.3, 

XTFRAC=1 .5,PITCH=  1 . 1 976,DSRA=.4920,DSLE=.0008,DSTE=.00 1 , 

NLE=48,NTE=36,WAKEP=1.,OMEGR=1.0,OMEGS=1.0,OMEGP=1.0,OMEGQ=1.0, 

&END 

&GRID3  AIRFX= 


127.0533, 

127.0025, 

126.9898, 

126.9467,  126.8984,  126.8400, 

126.7714, 

126.6927, 

126.5961, 

126.4768,  126.3752,  126.2990, 

126.2228, 

126.1466, 

125.9942, 

125.7656,  125.3084,  125.0798, 

124.9274, 

124.6988, 

124.5464, 

124.3940,  124.2416,  124.0892, 

124.0130, 

123.9368, 

123.8606, 

123.7844,  123.7082,  123.6320, 

123.5558, 

123.4796, 

123.3906, 

123.0276,  122.6646,  122.3016, 

121.9386, 

121.5757, 

121.2127, 

120.6435,  120.0744,  119.5052, 

118.9361, 

118.3669, 

117.7978, 

117.0168,  116.2359,  115.4549, 

114.6740, 

113.8930, 

113.1121, 

112.1089,  111.1057,  110.1025, 

109.0993, 

108.0961, 

107.0929, 

105.8614,  104.6300,  103.3985, 

102.1671, 

100.9356, 

99.7041, 

98.4585,  97.2128,  95.9671, 

94.7214, 

93.4757, 

92.2300, 

90.9738,  89.7175, 

88.4612, 

87.2049, 

85.9486, 

84.6923, 

83.4351,  82.1778, 

80.9205, 

79.6632, 

78.4059, 

77.1486, 

75.8955,  74.6423, 

73.3891, 

72.1360, 

70.8828, 

69.6296, 

68.3811,  67.1327, 

65.8842, 

64.6357, 

63.3872, 

62.1387, 

60.8997,  59.6607, 

58.4217, 

57.1827, 

55.9437, 

54.7047, 

53.4707,  52.2368, 

51.0028, 

49.7689, 

48.5350, 

47.3010, 

46.0678,  44.8346, 

43.6014, 

42.3682, 

41.1350, 

39.9019, 

38.6769,  37.4520, 

36.2270, 

35.0021, 

33.7771, 

32.5522, 

31.3293,  30.1064, 

28.8835, 

27.6606, 

26.4377, 

25.2148, 

24.0073,  22.7998, 

21.5923, 

20.3848, 

19.1774, 

17.9699, 

16.9731,  15.9764, 

14.9797, 

13.9829, 

12.9862, 

11.9894, 

11.2059,  10.4225, 

9.6390, 

8.8555, 

8.0720, 

7.2885,  t 

i.7058,  6.1230, 

5.5402, 

4.9574, 

4.3747, 

3.7919,  3 

1.4079,  3.0239, 

2.6399, 

2.2560, 

1.8720, 

1.4880,  1 

.4199,  1.3437, 

1.2675, 

1.1913, 

1.1151, 

1.0389,  0.9627,  0.8865, 

0.8103, 

0.7341, 

0.6579, 

0.5817,  0.5055,  0.4293, 

0.3531, 

0.2921, 

0.1956, 

0.1270,  0.0762,  0.0381, 

0.0102, 

-0.0051, 

-0.0127, 

-0.0152,  -0.0152,  -0.0127, 

-0.0051, 

0.0076, 

0.0279, 

0.0559,  C 

1.0889,  0.1270, 

0.1702, 

0.2210, 

0.2794, 

0.3454,  C 

1.4221,  0.7466, 

1.0712, 

1.3957, 

1.7202, 

2.0448,  2 

:.3693,  2.8953, 

3.4214, 
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3.9474, 

4.4734, 

4.9995, 

5.5255,  6 

.2765,  7.1 

0276, 

7.7786, 

8.5297, 

9.2807,  10.0317,  11.0113,  11.9910, 

12.9706, 

13.9503, 

14.9299, 

15.9096, 

17.1240, 

18.3385, 

19.5530, 

20.7675, 

21.9820, 

23.1965, 

24.4399, 

25.6833, 

26.9267, 

28.1700, 

29.4134, 

30.6568, 

31.9316, 

33.2064, 

34.4812, 

35.7560, 

37.0308, 

38.3056, 

39.6119, 

40.9182, 

42.2245, 

43.5307, 

44.8370, 

46.1433, 

47.4763, 

48.8092, 

50.1422, 

51.4752, 

52.8081, 

54.1411, 

55.4884, 

56.8358, 

58.1831, 

59.5304, 

60.8777, 

62.2250, 

63.5654, 

64.9057, 

66.2461, 

67.5864, 

68.9268, 

70.2671, 

71.5759, 

72.8847, 

74.1935, 

75.5023, 

76.8111, 

78.1199, 

79.3932, 

80.6665, 

81.9398, 

83.2130, 

84.4863, 

85.7596, 

87.0079, 

88.2561, 

89.5044, 

90.7527, 

92.0010, 

93.2493, 

94.4762, 

95.7031, 

96.9301, 

98.1570, 

99.3840, 

100.6109, 

101.8351, 

103.0593, 

104.2835, 

105.5077, 

106.7319, 

107.9562, 

108.9738, 

,  109.9914, 

111.0091, 

112.0267, 

113.0444, 

114.0620, 

114.8767, 

115.6915, 

116.5062, 

117.3210, 

118.1357, 

118.9505, 

119.5655, 

120.1806, 

120.7957, 

121.4108, 

122.0259, 

122.6410, 

123.0433, 

123.4456, 

123.8480, 

124.2503, 

124.6526, 

125.0549, 

125.2271, 

125.3795, 

125.5217, 

125.6487, 

125.7681, 

125.8799, 

125.9815, 

126.0805, 

126.1720, 

126.2583, 

126.3396, 

126.4133, 

126.4844, 

126.5504, 

126.6114, 

126.6673, 

126.7231, 

126.7714, 

126.8197, 

126.8628, 

126.9035, 

126.9390, 

126.9721, 

127.0025, 

127.0279, 

127.0508, 

127.0711, 

127.0864, 

127.1118, 

127.1245, 

127.1118, 

127.0787, 

127.0533, 

AIRFY= 


1.0617, 

0.9855, 

0.9093, 

0.8331, 

0.7569, 

0.6807, 

0.6045, 

0.5283, 

0.4521, 

0.3759, 

0.3200, 

0.2819, 

0.2489, 

0.2210, 

0.1702, 

0.1194, 

0.0838, 

0.0914, 

0.1067, 

0.1448, 

0.1803, 

0.2210, 

0.2667, 

0.3200, 

0.3505, 

0.3810, 

0.4140, 

0.4470, 

0.4826, 

0.5232, 

0.5613, 

0.6045, 

0.6577, 

0.8979, 

1.1299, 

1.3544, 

1.5718, 

1.7825, 

1.9871, 

2.2968, 

2.5944, 

2.8820, 

3.1614, 

3.4345, 

3.7031, 

4.0671, 

4.4257, 

4.7782, 

5.1240, 

5.4623, 

5.7926, 

6.2042, 

6.6019, 

6.9862, 

7.3577, 

7.7168, 

8.0641, 

8.4753, 

8.8712, 

9.2538, 

9.6248, 

9.9861, 

10.3395, 

10.6901, 

11.0331, 

11.3669, 

11.6900, 

12.0012, 

12.2988, 

12.5843, 

12.8558, 

13.1147, 

13.3621, 

13.5993, 

13.8275, 

14.0479, 

14.2594, 

14.4611, 

14.6519, 

14.8307, 

14.9963, 

15.1474, 

15.2840, 

15.4055, 

15.5113, 

15.6010, 

15.6742, 

15.7302, 

15.7704, 

15.7957, 

15.8072, 

15.8060, 

15.7932, 

15.7698, 

15.7350, 

15.6877, 

15.6268, 

15.5514, 

15.4603, 

15.3532, 

15.2293, 

15.0884, 

14.9302, 

14.7545, 

14.5611, 

14.3502, 

14.1228, 

13.8804, 

13.6245, 

13.3564, 

13.0778, 

12.7917, 

12.4966, 

12.1926, 

11.8797, 

11.5580, 

11.2275, 

10.8891, 

10.5426, 

10.1886, 

9.8278, 

9.4608, 

9.0881, 

8.7149, 

8.3365, 

7.9522, 

7.5617, 

7.1646, 

6.7604, 

6.4211, 

6.0769, 

5.7279, 

5.3741, 

5.0158, 

4.6531, 

4.3647, 

4.0722, 

3.7742, 

3.4694, 

3.1562, 

2.8332, 

2.5860, 

2.3330, 

2.0747, 

1.8115, 

1.5438, 

1.2721, 

1.0910, 

0.9086, 

0.7247, 
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0.5397, 

0.3536, 

0.1666, 

0.0457, 

0.0279, 

0.0127, 

-0.0025, 

0.0025, 

0.0178, 

0.1346, 

0.2108, 

0.2870, 

0.5918, 

0.6680, 

0.7391, 

0.9728, 

1.0490, 

1.1252, 

1.4300, 

1.5062, 

1.5824, 

2.5977, 

2.8950, 

3.1850, 

4.7579, 

5.1599, 

5.5503, 

7.4587, 

7.9337, 

8.3922, 

10.4250, 

10.9137, 

11.3861, 

13.4531, 

13.9565, 

14.4440, 

16.2726, 

16.6967, 

17.1081, 

18.6653, 

19.0255, 

19.3718, 

20.6277, 

20.9017, 

21.1577, 

22.0006, 

22.1587, 

22.2936, 

22.5896, 

22.6023, 

22.5911, 

22.3138, 

22.1843, 

22.0292, 

21.1691, 

20.8976, 

20.6055, 

19.2894, 

18.9281, 

18.5564, 

17.0114, 

16.6061, 

16.1909, 

14.4552, 

14.0081, 

13.5599, 

11.7911, 

11.3519, 

10.9113, 

9.3460, 

8.9668, 

8.5851, 

7.2734, 

6.9654, 

6.6594, 

5.6775, 

5.4485, 

5.2154, 

4.4871, 

4.3159, 

4.1398, 

3.7287, 

3.6525, 

3.5763, 

3.2715, 

3.1953, 

3.1191, 

2.8143, 

2.7381, 

2.6619, 

2.3571, 

2.2809, 

2.2047, 

1.8999, 

1.8237, 

1.6713, 

1.0617, 

&END 


0.1321, 

0.0991, 

0.0711, 

0.0000, 

-0.0025, 

-0.0051, 

0.0381, 

0.0660, 

0.0991, 

0.3632, 

0.4394, 

0.5156, 

0.7442, 

0.8204, 

0.8966, 

1.2014, 

1.2776, 

1.3538, 

1.6596, 

1.9804, 

2.2930, 

3.4682, 

3.9135, 

4.3429, 

5.9304, 

6.4574, 

6.9667, 

8.8345, 

9.3880, 

9.9173, 

11.8448, 

12.3975, 

12.9336, 

14.9158, 

15.3830, 

15.8350, 

17.5077, 

17.9059, 

18.2919, 

19.7037, 

20.0282, 

20.3364, 

21.3953, 

21.6180, 

21.8200, 

22.4043, 

22.4910, 

22.5527, 

22.5566, 

22.4993, 

22.4186, 

21.8481, 

21.6457, 

21.4189, 

20.2942, 

19.9741, 

19.6386, 

18.1759, 

17.7956, 

17.4077, 

15.7652, 

15.3362, 

14.8987, 

13.1131, 

12.6704, 

12.2301, 

10.4680, 

10.0967, 

9.7226, 

8.2011, 

7.8922, 

7.5826, 

6.3565, 

6.1300, 

5.9042, 

4.9770, 

4.8173, 

4.6541, 

3.9587, 

3.8811, 

3.8049, 

3.5001, 

3.4239, 

3.3477, 

3.0429, 

2.9667, 

2.8905, 

2.5857, 

2.5095, 

2.4333, 

2.1285, 

2.0523, 

1.9761, 

1.4808, 

1.2903, 

1.1379, 
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“RVCQ3D”  INPUT 


'GELDER  CONTROLLED  DIFFUSION  BLADE' 

&NLI  M=340,N=49,MTL=80,MIL=143  &end 

&NL2  NSTG=4,IVDT=LIRS=1,EPI=.30,EPJ=.40,CFL=4.0  AVISC2=0.,  EPS=0.4, 
AVISC4=1.  &END 

&NL3  IBCIN=l,IBCEX=l,ITMAX=5000,IRESTI=0,IRESTO=l,IRES=l,ICRNT=I0000, 
IXRM=0&END 

&NL4  AMLE=0.22,ALLE=36.5,BETE=0,PRAT=O.9767,P0IN=I .0,T0IN=1 .0,G=I .40, 
&END 

&NL5  ILT=4, JEDGE=30,RENR=  1 .0e6,PRNR=.70,TW=0.0, V1SPWR=.666666, 
CMUTM=14.0.ITUR=2  &END 

&NL6  OMEGA=0.0,NBLADE=1,NMN=2  &end 

&NL7  TINTENS=0.02  &end 
-.3000  3.000 
1.0000  1.000 
1.0000  .9709 
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“PCP”  CODE 


c - c 

C  PROGRAM  TO  READ  THE  OUTPUT  FROM  RVCQ3D.F  AND  GRAPE.F  C 
C  AND  GENERATE  A  DATA  FILE  FOR  A  P/PO  VS  CHORD  PLOT  C 
C  AROUND  THE  AIRFOIL.  AUG.22,  1991  C 

C - C 

c  IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  Q(400,400,4),X(400, 400),  Y(400, 400) 

DIMENSION  U(400,400),V(400,400),P(400,400) 

C 

C  CALCULATE  THE  CRITICAL  VELOCITY 
cepe=1.00 
rgas=1.0 
pi=1.0 
ti=L0 

PRAT=0.685 

rhoi=pi/(rgas*ti) 

ceve=cepe-rgas 

g=1.40 

gP=g+l 

cstar=sqrt(2  *g’''pi/(gp*rhoi)) 

C 

ISTART=80 

IFINIT=261 

C 

READ(7,*)NI,NJ 
print  *,ni,nj 

READ(7,*)((X(I,J),I=1,NI),J=1,NJ), 

&  ((Y(I,J),I=1.NI),J=1,NJ) 

READ(7,*)MTL,MIL 

C 

READ(3,*)NI,NJ 

READ(3,*)FSMACH,ALF,RE,TIME 

READ(3,*)(((Q(I,J,K),I=1,NI),J=1,NJ),K=1,4) 

C 

GAMMA=1.4 

C 

DO  1234  I=1,NI 
D0  4321  J=1,NJ 
U(I,J)=Q(1,J,2)/Q(I,J,1) 

V(I,J)=Q(I,J,3)/Q(I,J,1) 

P(I,J)=(GAMMA-1)*(Q(I,J,4)-.5*Q(I,J,1)*(U(I,J)**2+V(I,J)**2)) 

&  ♦RHOI*CSTAR**2 
4321  CONTINUE 
1234  CONTINUE 
C 

C  DETERMINE  THE  XMIN  AND  XMAX  GRID  POINT  POSITION 
XMIN=0.0 
XMAX=0.0 
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DO  2  I=ISTART,IFINIT 
IF(X(1, 1  ).LE.XMIN)THEN 
IMIN=I 
PRINT  *,IMIN 
XMIN=X(I,1) 

END  IF 

IF(X(1, 1).GE.XMAX)THEN 
IMAX=I 

c  PRINT  *,IMAX 
XMAX=X(I,I) 

END  IF 
2  CONTINUE 

FREE  STREAM  STATIC  PRESSURE  USES  THE  XMIN  GRID  POINT 
P0=P(IMIN,NJ)+.5*Q(IMIN,NJ,1)*(U(IMIN,NJ)**2+V(IMIN,NJ)**2) 

&  *RHOI*CSTAR**2 

i=imin 
j=nj 

lO=l/rgas*((g-I)*q(iJ,4)/q(iJ,lHg-l)**2/(2*g)* 

&  (q(ij,2)**2+q(ij,3)**2)/q(iJ,l)*-^2) 

tO=lO*cstar**2 

print  *,pO,lO,p(imin,nj),q(imin,nj,l),u(iinin,nj),v(imin,nj) 

CALCULATE  THE  MASS  FLOW  RATE  FROM  THE  INPUT  DATA 
NORMALIZED  WITH  THE  INLET  AREA 
dmass=pi/(rgas*ti)*(l/prat)**(-l/g)*sqrt(g*(g-l)*rgas*ti/2* 

&  (l-I/(l/prat)**((g-l)/g))) 

print  *,'cal.  inlet  mass  flow  rate  =',dmass,'*  area' 


dmasse=q(imin,nj,l)’'‘RHOI*u(imin,nj) 

print  *,'comp.  inlet  mass  flow  rate  /  area  =',dmasse 

PRINT  OUT  DOWN  STREAM  CONDITIONS 
pl=p(l,l)+.5*q(Ll,l)*(u(Ll)*'*2+v(Ul)**2)*RHOI*CSTAR**2 
tl=l/rgas*((g-l)*q(LL4)/q(Ll,l)-(g-l)**2/(2*g)* 

&  (q(l,l,2)**2+q(l,l,3)**2)/q(LLl)**2) 
tl=tl*CSTAR**2 
print  ♦,'downstream  condition' 
print  *,pLtLp(Ll),q(Ll,l),u(l,l),v(l,l) 

calculate  the  mass  flow  rate  at  exit 
smass=0.0 
ra=0.0 
ba=0.0 

do  21  j=l,nj-l 

smass=smass+(q(  1  j,  1  )+q(  1 J+ 1 , 1  ))*(u(  1  J)+u(  1 ,  J+ 1  ))* 

&  RHOI*CSTAR*(abs(y(  1  j+ l)-y(  1  j)))*0.25 

ra=ra+(q(lj,l)+q(lj+l,l))*(abs(y(lj+l)-y(lj)))*0.5* 

&  RHOI 

ba=ba+abs(y(  1  j+ 1  )-y(  1  j)) 
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21  continue 

do  22  j=l,nj-l 

smass=smass+(q(ni  j,  1  )+q(ni  j+ 1 , 1  ))*(u(ni  j)+u(ni,  J+ 1  ))♦ 

&  RHOI*CSTAR*(abs(y(ni  j+ 1  )-y(ni  J)))*0.25 

ra=ra+(q(n  1  j,  1  )+q(ni  J+ 1 , 1  ))*(abs(y(ni  J+ l)-y(ni  J)))*0.5 
&  *RHOI 

ba=ba+abs(y(nij+l)-y(nij)) 

22  continue 

average  velocity 
va=smass/ra 
dr=ra^a 

print  *>61.  downsteam,vd  =',va 
print  *,'ave.  density,dr  -  ,dr 

CALCULATE  THE  INPUT  INLET  VELOCITY 
NORMALIZED  BY  CSTAR 
cO=sqrt(g*rgas*ti) 

vin=sqrt(g*(g-l)*rgas*ti/2*(l-(l/prat)**((l-g)/g)))/c0 
print  *,'input  inlet  velocity  /  cO  =’,vin 

CALCULATE  THE  CORRECTED  DOWNSTREAM  VELOCITY 
do23  j=nj,l,-l 
vac=u(nij)*vin/va 
alpha=atan(abs(v(nij)/u(nij))) 
vtc=vac*lan(alpha) 
vc=sqrt(vac**2+vtc**2) 
yy=abs(y(ni,nj)-y(Lnj)) 
write(24,*)abs(y(nij)-y(ni,nj))/yy,vc 

23  continue 

do  25  j=2,nj 
vac=u(lj)*vin/va 
alpha=atan(abs(v(  1  j)/u(  1  j))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac*  *2+vtc*  *2) 
write(24,*)abs(y(Ij)-y(ni,nj))/yy,vc 
25  continue 


CHORD=ABS(X(IM  AX,  1  )-X(IMIN,  1 )) 

C 

DO  1  I=ISTART,IFINIT 
DIST=ABS(X(I,  I  )-X(IMIN,  1 )) 

XS=DIST/CHORD 

CP=I.O+(-PO+P(I,l))/(.5*Q(IMIN,NJ,l)*(U(IMIN,NJ)**2 
&  +V(IMIN,NJ)**2)*RHOI*CSTAR**2) 

vvcr=sqrt(u(i,  l)**2+v(i,l)**2) 

C  THE  SURFACE  POINT  P(1, 1 )  ALWAYS  EQUAL  0  ? 

PT=P(I,1)/P0 

c  WRITE(63,*)XS,WCR 
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c  WRITE(64,*)XS,PT 
WRITE(65,*)XS,CP 
1  CONTINUE 
C 

STOP 

END 
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